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Based  on  concepts  drawn  from  the  ekpyrotic  scenario  and  M-theory,  we  elab¬ 
orate  onr  recent  proposal  of  a  cyclic  model  of  the  Universe.  In  this  model, 
the  Universe  undergoes  an  endless  sequence  of  cosmic  epochs  which  begin 
with  the  Universe  expanding  from  a  ‘big  bang’  and  end  with  the  Universe 
contracting  to  a  ‘big  crunch.’  Matching  from  ‘big  crunch’  to  ‘big  bang’  is  per¬ 
formed  according  to  the  prescription  recently  proposed  with  Khoury,  Ovrut 
and  Seiberg.  The  expansion  part  of  the  cycle  includes  a  period  of  radiation 
and  matter  domination  followed  by  an  extended  period  of  cosmic  accelera¬ 
tion  at  low  energies.  The  cosmic  acceleration  is  crucial  in  establishing  the 
flat  and  vacuous  initial  conditions  required  for  ekpyrosis  and  for  removing 
the  entropy,  black  holes,  and  other  debris  produced  in  the  preceding  cycle. 
By  restoring  the  Universe  to  the  same  vacuum  state  before  each  big  crunch, 
the  acceleration  insures  that  the  cycle  can  repeat  and  that  the  cyclic  solution 
is  an  attractor. 
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1  Introduction 


In  a  recent  paper, 0  we  introduced  the  possibility  of  a  cyclic  Universe,  a  cos¬ 
mology  in  which  the  Universe  undergoes  a  periodic  sequence  of  expansion  and 
contraction.  Each  cycle  begins  with  a  “big  bang”  and  ends  in  a  “big  crunch,” 
only  to  emerge  in  a  big  bang  once  again.  The  expansion  phase  of  each  cycle 
includes  a  period  of  radiation-,  matter-,  and  quintessence-domination,  the 
last  phase  of  which  corresponds  to  the  current  epoch  of  cosmic  acceleration. 
The  accelerated  expansion  phase  dilutes  by  an  exponential  factor  the  en¬ 
tropy  and  the  density  of  black  holes  and  any  other  debris  produced  since  the 
preceding  big  bang.  The  acceleration  ultimately  ends,  and  it  is  followed  by 
a  period  of  decelerating  expansion  and  then  contraction.  At  the  transition 
from  big  crunch  to  big  bang,  matter  and  radiation  are  created,  restoring  the 
Universe  to  the  high  density  required  for  a  new  big  bang  phase. 

Historically,  cyclic  models  have  been  considered  attractive  because  they 
avoid  the  issue  of  initial  conditions^  Examples  can  be  found  in  mythologies 
and  philosophies  dating  back  to  the  beginning  of  recorded  history.  Since 
the  introduction  of  general  relativity,  though,  various  problems  with  the 
cyclic  concept  have  emerged.  In  the  1930 ’s,  Richard  Tolmani  discussed  cyclic 
models  consisting  of  a  closed  Universe  with  zero  cosmological  constant.  He 
pointed  out  that  entropy  generated  in  one  cycle  would  add  to  the  entropy 
created  in  the  next.  Consequently,  the  maximal  size  of  the  Universe,  and 
the  duration  of  a  cycle,  increase  from  bounce  to  bounce.  Extrapolating 
backwards,  the  duration  of  the  bounce  converges  to  zero  in  a  finite  time. 
Consequently,  the  problem  of  initial  conditions  remains.  In  the  1960’s,  the 
singularity  theorems  of  Hawking  and  Penrose  showed  that  a  big  crunch  neces¬ 
sarily  leads  to  a  cosmic  singularity  where  general  relativity  becomes  invalid. 
Without  a  theory  to  replace  general  relativity  in  hand,  considerations  of 
whether  time  and  space  could  exist  before  the  big  bang  were  discouraged. 
“Big  bang”  became  synonymous  with  the  beginning  of  space-time.  However, 
there  is  nothing  in  the  Hawking-Penrose  singularity  theorems  to  suggest  that 
cyclic  behavior  is  forbidden  in  an  improved  theory  of  gravity,  such  as  string 
theory  and  M  theory,  and  some  people  have  continued  to  speculate  on  this 
possibility.!’ I  In  the  1990’s,  observations  showed  that  the  matter  density  is 
significantly  less  than  the  critical  density  and  that  the  scale  factor  of  the 
Universe  is  accelerating.!  Tolman’s  cyclic  model  based  on  a  closed  Universe 
is  therefore  observationally  ruled  out. 
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Curiously,  the  same  observations  that  eliminate  Tolman’s  cyclic  model 
fit  perfectly  the  novel  kind  of  cyclic  model  proposed  here.  In  our  proposal, 
the  Universe  is  flat,  rather  than  closed.  The  transition  from  expansion  to 
contraction  is  caused  by  introducing  negative  potential  energy,  rather  than 
spatial  curvature.  Furthermore,  the  cyclic  behavior  depends  in  an  essen¬ 
tial  way  on  having  a  period  of  accelerated  expansion  after  the  radiation  and 
matter-dominated  phases.  During  the  accelerated  expansion  phase,  the  Uni¬ 
verse  approaches  a  nearly  vacuous  state,  restoring  very  nearly  identical  local 
conditions  as  existed  in  the  previous  cycle  prior  to  the  contraction  phase. 
Globally,  the  total  entropy  in  the  Universe  grows  from  cycle  to  cycle,  as  Tol- 
man  suggested.  However,  the  entropy  density,  which  is  all  any  real  observer 
would  actually  see,  has  perfect  cyclic  behavior  with  entropy  density  being 
created  at  each  bounce,  and  subsequently  being  diluted  to  negligible  levels 
before  the  next  bounce. 

The  linchpin  of  the  cyclic  picture  is  safe  passage  through  the  cosmic  sin¬ 
gularity,  the  transition  from  the  big  crunch  to  big  bang.  In  recent  work  with 
J.  Khoury,  B.  Ovrut  and  N.  Seiberg,  we  have  proposed  that  a  smooth  tran¬ 
sition  is  possible  in  string  theory.Q’U  In  ordinary  4d  general  relativity,  the  big 
crunch  is  interpreted  as  the  collapse  and  disappearance  of  four- dimensional 
space-time.  Densities  and  curvatures  diverge  and  and  there  is  no  sign  that 
a  transition  is  possible.  But  in  the  theory  considered  here,  what  appears  to 
be  a  big  crunch  in  the  4d  effective  theory  actually  corresponds  to  the  mo¬ 
mentary  collapse  of  an  additional  fifth  dimension.  As  far  as  matter  which 
couples  to  the  higher  dimensional  metric  is  concerned,  the  three  large  spa¬ 
tial  dimensions  remain  large  and  time  continues  smoothly.  The  temperature 
and  density  are  finite  as  one  approaches  the  crunch,  and,  furthermore,  the 
geometry  is  flat  just  before  and  just  after  the  bounce.  In  short,  there  is 
nothing  to  suggest  that  time  comes  to  an  end  when  the  fifth  spatial  dimen¬ 
sion  collapses.  Quite  the  contrary,  the  most  natural  possibility  is  that  time 
continues  smoothly.  Efforts  are  currently  underway  to  establish  this  conclu¬ 
sion  rigorously  in  string  theory.0  The  cyclic  scenario  considered  here  exploits 
this  concept  and  is  absolutely  dependent  on  its  validity.  In  the  absence  of 
a  detailed  theory  of  the  transition  from  big  crunch  to  big  bang,  we  will  pa¬ 
rameterize  the  bounce  in  terms  of  simple  matching  conditions  incorporating 
energy  and  momentum  conservation. 

The  appeal  of  a  cyclic  model  is  that  it  provides  a  description  of  the  history 
of  the  Universe  which  applies  arbitrarily  far  back  into  our  past.  The  model 
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presented  here  suggests  novel  answers  to  some  of  the  most  challenging  issues 
in  cosmology:  How  old  is  the  Universe  -  finite  or  infinite?  How  large  is 
it?  What  was  the  physical  cause  of  its  homogeneity,  isotropy  and  flatness? 
What  was  the  origin  of  the  energy  density  inhomogeneities  that  seeded  cosmic 
structure  formation  and  are  visible  on  the  cosmic  microwave  sky?  What  is 
the  resolution  of  the  cosmic  singularity  puzzle?  Was  there  time,  and  an 
arrow  of  time,  before  the  big  bang?  In  addition,  our  scenario  has  a  number 
of  surprising  implications  for  other  major  puzzles  such  as  the  value  of  the 
cosmological  constant,  the  relative  densities  of  different  forms  of  matter,  and 
even  for  supersymmetry  breaking. 

The  cyclic  model  rests  heavily  on  ideas  developed  as  part  of  the  recently 
proposed  “ekpyrotic  Universe. ”01  The  basic  physical  notion  is  that  the  col¬ 
lision  between  two  brane  worlds  approaching  one  another  along  an  extra 
dimension  would  have  literally  generated  a  hot  big  bang.  Although  the  orig¬ 
inal  ekpyrosis  paper  focused  on  collisions  between  bulk  branes  and  boundary 
branes, 0  here  the  more  relevant  example  is  where  the  boundary  branes  col¬ 
lide,  the  extra  dimension  disappears  momentarily  and  the  branes  then  bounce 
apart!  The  ekpyrotic  scenario  introduced  several  important  concepts  that 
serve  as  building  blocks  for  the  cyclic  scenario: 

•  boundary  branes  approaching  one  another  (beginning  from  rest)  corre¬ 
sponds  to  contraction  in  the  effective  4d  theoretic  description;! 

•  contraction  produces  a  blue  shift  effect  that  converts  gravitational  en¬ 
ergy  into  brane  kinetic  energy;! 

•  collision  converts  some  fraction  of  brane  kinetic  energy  into  matter  and 
radiation  that  can  fuel  the  big  bang; 

•  the  collision  and  bouncing  apart  of  boundary  branes  corresponds  to  the 
transition  from  a  big  crunch  to  a  big  bangJa 

A  key  element  is  added  to  obtain  a  cyclic  Universe.  The  ekpyrotic  scenario 
assumes  that  there  is  only  one  collision  after  which  the  interbrane  potential 
becomes  zero  (perhaps  due  to  changes  in  the  gauge  degrees  of  freedom  on  the 
branes  that  zero  out  the  force).  The  cyclic  model  assumes  instead  that  the 
interbrane  potential  is  the  same  before  and  after  collision.  After  the  branes 
bounce  and  fly  apart,  the  interbrane  potential  ultimately  causes  them  to  draw 
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together  and  collide  again.  To  ensure  cyclic  behavior,  we  will  show  that  the 
potential  must  vary  from  negative  to  positive  values.  (In  the  ekpyrotic  ex¬ 
amples,  the  potentials  are  zero  or  negative  for  all  inter-brane  separations.) 
We  propose  that,  at  distances  corresponding  to  the  present-day  separation 
between  the  branes,  the  inter-brane  potential  energy  density  should  be  posi¬ 
tive  and  correspond  to  the  currently  observed  dark  energy,  providing  roughly 
70%  of  the  critical  density  today.  That  is,  the  dark  energy  that  is  causing 
the  cosmic  acceleration  of  the  Universe  today  is,  in  this  scenario,  inter-brane 
potential  energy.  The  dark  energy  and  its  associated  cosmic  acceleration  play 
an  essential  role  in  restoring  the  Universe  to  a  nearly  vacuous  state  thereby 
allowing  the  cyclic  solution  to  become  an  attractor.  As  the  brane  separation 
decreases,  the  interbrane  potential  becomes  negative,  as  in  the  ekpyrotic  sce¬ 
nario.  As  the  branes  approach  one  another,  the  scale  factor  of  the  Universe, 
in  the  conventional  Einstein  description,  changes  from  expansion  to  contrac¬ 
tion.  When  the  branes  collide  and  bounce,  matter  and  radiation  are  produced 
and  there  is  a  second  reversal  transforming  contraction  to  expansion  so  a  new 
cycle  can  begin. 

The  central  element  in  the  cyclic  scenario  is  a  four  dimensional  scalar 
held  0,  parameterizing  the  inter-brane  distance  or  equivalently  the  size  of 
the  fifth  dimension.  The  brane  separation  goes  to  zero  as  0  tends  to  —  oo, 
and  the  maximum  brane  separation  is  attained  at  some  finite  value  0max. 
For  the  most  part  our  discussion  will  be  framed  completely  within  the  four 
dimensional  effective  theory  of  gravity  and  matter  coupled  to  the  scalar  held 
0.  This  description  is  universal  in  the  sense  that  many  higher  dimensional 
brane  models  converge  to  the  same  four  dimensional  effective  description  in 
the  limit  of  small  brane  separation.  We  shall  not  need  to  tie  ourselves  to 
a  particular  realization  of  the  brane  world  idea,  such  as  heterotic  M  theory 
for  the  purposes  of  this  discussion,  although  such  an  underlying  description 
is  certainly  required,  both  for  actually  deriving  the  scalar  potential  we  shall 
simply  postulate  and  for  the  ultimate  quantum  consistency  of  the  theory. 
The  extra  dimensional,  and  string  theoretic  interpretation,  is  also  crucial 
at  the  brane  collision,  where  the  effective  four  dimensional  Einstein- frame 
description  is  singular  and  at  which  point  we  postulate  a  big  crunch-big 
bang  transition  as  outlined  in  Ref.  Again,  for  the  present  discussion  we 
simply  parameterize  the  outcome  of  this  transition  in  terms  of  the  density 
of  radiation  produced  on  the  branes,  and  the  change  in  the  kinetic  energy  of 
the  scalar  held,  corresponding  to  a  change  in  the  contraction/expansion  rate 
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of  the  fifth  dimension. 

The  scalar  held  0  plays  a  crucial  role  in  the  cyclic  scenario,  in  regularizing 
the  Einstein-frame  singularity.  Matter  and  radiation  on  the  brane  couple  to 
the  Einstein  frame  scale  factor  a  times  a  function  0(0)  with  exponential 
behavior  as  0  — >  —  oo,  such  that  the  product  is  generically  finite  at  the  brane 
collision,  even  though  a  =  0  and  0  =  —  oo  there.  For  finite  0,  the  coupling 
of  the  matter  and  radiation  to  0  is  more  model-dependent.  Models  in  which 
0  is  massless  at  the  current  epoch,  such  as  we  describe  in  this  paper,  face 
a  strong  constraint  due  to  the  fact  that  0  can  mediate  a  'fifth  force’,  which 
is  in  general  composition  dependent  and  violates  the  equivalence  principle. 
Again,  without  tying  ourselves  to  a  particular  brane  world  scenario  we  shall 
consider  models  in  which  the  coupling  function  0(0)  tends  to  a  constant  at 
current  values  of  0  (large  brane  separations),  and  the  corresponding  fifth 
force  is  weak.  An  example  of  such  a  model  is  the  Randall-Sundrum  model 
with  the  non-relativistic  matter  we  are  made  of  localized  on  the  positive 
tension  brane  (see  e.g.  Ref.  [lO]  for  a  recent  discussion).  In  models  where 


0(0)  does  not  tend  to  a  constant  at  current  values  of  0,  one  must  invoke 
some  physical  mechanism  to  give  the  0  held  a  small  mass  so  that  the  fifth 
force  is  only  short-ranged.  This  modification  still  allows  for  cyclic  behavior, 
with  an  epoch  of  false  vacuum  domination  followed  by  tunneling.0 

The  outline  of  this  paper  is  as  follows.  In  Section  ^],  we  describe  the  req¬ 
uisite  properties  of  the  scalar  held  (inter-brane)  potential  and  present  a  brief 
tour  through  one  complete  cosmic  cycle.  In  subsequent  sections,  we  focus  in 
technical  detail  on  various  stages  of  the  cycle:  the  bounce  (Section  [5|),  pass¬ 
ing  through  the  potential  well  after  the  big  bang  (Section  fD,  the  radiation-, 
matter-  and  quintessence-  dominated  epochs  (Section  ||),  the  onset  of  the 
contraction  phase  and  the  generation  of  density  perturbations  (Section  ||). 
In  Section  [!],  we  show  that  the  cyclic  solution  is  a  stable  attractor  solution 
under  classical  and  quantum  huctuations.  In  Section  |],  we  discuss  the  im¬ 
plications  for  the  fundamental  questions  of  cosmology  introduced  above. 


2  A  Brief  Tour  of  the  Cyclic  Universe 


The  various  stages  of  a  cyclic  model  can  be  characterized  in  terms  of  a  scalar 
held  0  which  moves  back  and  forth  in  an  effective  potential  V(0).  In  Sec¬ 
tion  p~T|.  we  discuss  the  basic  properties  that  V(0)  must  have  in  order  to 
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allow  cyclic  solutions. 

The  stages  of  expansion  and  contraction  can  be  described  from  two  dif¬ 
ferent  points  of  view.  First,  one  can  choose  fields  and  coordinates  so  that 
the  full  extra-dimensional  theory  is  reduced  to  an  effective  four- dimensional 
theory  with  a  conventional  Einstein  action.  The  key  parameters  are  the  scale 
factor  a  and  the  modulus  scalar  field  0  that  determines  the  distance  between 
branes.  In  this  picture,  the  terms  “big  bang”  and  “big  crunch”  seem  well- 
merited.  The  scale  factor  collapses  to  zero  at  the  big  crunch,  bounces,  and 
grows  again  after  the  big  bang.  However,  what  is  novel  is  the  presence  of 
the  scalar  field  0  which  is  runs  to  — oo  at  the  bounce  with  diverging  kinetic 
energy.  The  scalar  field  acts  as  a  fifth  force,  modifying  in  an  essential  way 
the  behavior  of  matter  and  energy  at  the  big  crunch.  Namely,  the  tempera¬ 
ture  and  matter  density  remain  finite  at  the  bounce  because  the  usual  blue 
shift  effect  during  contraction  is  compensated  by  the  fifth  force  effect  due 
to  0.  The  arrangement  seems  rather  magical  if  one  is  unaware  that  the  4d 
theory  is  derived  from  a  higher  dimensional  picture  in  which  this  behavior 
has  a  clear  geometrical  interpretation.  Nevertheless,  for  most  of  this  paper 
we  shall  keep  to  the  four  dimensional  Einstein  description,  switching  to  the 
higher  dimensional  picture  only  when  necessary  to  understand  the  bounce, 
or  to  discuss  global  issues  where  matching  one  cycle  to  the  next  is  important. 
The  description  of  a  single  cycle  from  the  4d  effective  theory  point-of-view  is 
given  in  Section  [2.2[. 

The  same  evolution  appears  to  be  quite  different  to  observers  on  the  vis¬ 
ible  brane  who  detect  matter  and  radiation  confined  to  three  spatial  dimen¬ 
sions.  In  this  picture,  depending  on  the  details,  the  brane  is  either  always,  or 
nearly  always  expanding  except  for  tiny  jags  near  the  big  crunch/big  bang 
transition  when  it  contracts  by  a  modest  amount.  The  branes  stretch  at 
a  rate  that  depends  on  which  form  of  energy  dominates  the  energy  density 
of  the  Universe.  As  the  big  crunch  is  approached,  however,  the  expansion 
rate  changes  suddenly,  and  new  matter  and  radiation  is  created  (a  brane  has 
instantaneously  collided  into  the  visible  brane  and  bounced  from  it).  We 
describe  some  aspects  of  the  visible  brane  viewpoint  in  Section  |2.3|. 

This  picture  makes  it  clear  that  the  big  crunch  does  not  correspond  to 
the  disappearance  of  all  of  space  and  the  end  of  time  but,  rather,  to  the 
momentary  disappearance  of  a  fifth  dimension.  However,  the  behavior  of 
gravity  itself  appears  quite  wild  because  it  depends  on  the  full  bulk  space- 
time,  which  is  changing  rapidly.  One  way  of  describing  this  picture  is  that 
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one  has  mapped  the  conventional  big  bang  singularity  onto  the  mildest  form 
of  singularity  possible,  namely  the  disappearance  of  a  single  dimension  for 
an  instant  of  time.  Nevertheless  there  are  delicate  issues  involved,  as  are  dis¬ 
cussed  in  Ref.  [|,  such  as  the  fact  that  the  effective  four  dimensional  Planck 
mass  hits  zero  at  the  singularity,  so  that  gravitational  fluctuations  can  be¬ 
come  large.  There  are  suggestions  in  specific  calculation^!!  that  physical 
quantities  are  nevertheless  well  behaved  although  a  great  deal  more  remains 
to  be  done  to  make  the  picture  rigorous. 

2.1  The  Effective  Potential  for  a  Cyclic  Universe 

We  will  consider  in  this  paper  potentials  U(0)  of  the  form  shown  in  Figure  1, 
with  the  following  key  features: 

•  The  potential  tends  to  zero  rapidly  as  0  — »  —  oo.  One  natural  possibility 
for  the  extra  dimension  parameterized  by  0  is  the  eleventh  dimension  of  M 
theory.  In  this  case  the  string  coupling  constant  gs  oc  e7<y  with  some  positive 
constant  7,  and  gs  vanishes  as  0  — >  —00.  Non-perturbative  potentials  should 
vanish  faster  than  any  finite  power  of  gs,  i.e.,  faster  than  an  exponential  in 
4>. 

•  The  potential  is  negative  for  intermediate  0,  and  rises  with  a  region  of 
large  negative  curvature,  V" jV  >>  1  covering  a  range  of  0  of  order  unity 
in  Planck  mass  units.  This  region  is  required  for  the  production  of  scale 
invariant  density  perturbations,  as  proposed  in  Ref.  [7]  and  detailed  in  Ref.  [T2|. 
Attractive  exponential  potentials  of  this  type  could  be  produced,  for  example, 
by  the  virtual  exchange  of  massive  particles  between  the  boundary  branes. 

•  As  0  increases,  the  potential  rises  to  a  shallow  plateau,  with  V"  jV  <<  1 
and  a  positive  height  Vo  given  by  the  present  vacuum  energy  of  the  Universe 
as  inferred  from  cosmic  acceleration  and  other  astronomical  evidence.  The 
positive  energy  density  is  essential  for  having  a  cyclic  solution  since  it  pro¬ 
duces  a  period  of  cosmic  acceleration  that  restores  the  Universe  to  a  nearly 
vacuous  state  before  the  next  bounce.  The  discussion  here  can  be  extended 
to  potentials  of  a  more  general  form.  For  example,  it  is  not  essential  that 
the  positive  plateau  persist  to  arbitrarily  large  0  since  the  cyclic  solution 
only  explores  a  finite  range  of  0  >  0.  Provided  the  condition  V" /V  -C  1  is 
satisfied  over  that  range,  the  Universe  undergoes  cosmic  acceleration  when 
the  held  rolls  down  that  portion  of  the  potential.  However,  for  simplicity,  we 
will  consider  the  example  in  Figure  1. 


An  explicit  model  for  V (0)  which  is  convenient  for  analysis  is 


v(4>)  =  Vi(l  -  c-*)F(*), 


(1) 


where,  without  loss  of  generality,  we  have  shifted  0  so  that  the  zero  of  the 
potential  occurs  at  0  =  0.  The  function  U(0)  is  introduced  to  represent 
the  vanishing  of  non-perturbative  effects  described  above:  F(0)  turns  off  the 
potential  rapidly  as  0  goes  below  (pmin,  but  it  approaches  one  for  0  >  0min. 
For  example,  F(0)  might  be  proportional  to  e~lt9s  or  e”1//9s,  where  gs  oc  e7^ 
for  7  >  0.  The  constant  Vo  is  set  roughly  equal  to  the  vacuum  energy 
observed  in  today’s  Universe,  of  order  10“ 120  in  Planck  units.  We  do  not 
attempt  to  explain  this  number.  Various  suggestions  as  to  how  a  suitable 
small  positive  vacuum  energy  could  arise  have  been  made.citesw,bp  For  large 
c,  this  potential  has  V" /V  <C  1  for  0  ^  1  and  V" /V  1  for  0miri  < 
0  <  0.  These  two  regions  account  for  cosmic  acceleration  and  for  ekpyrotic 
production  of  density  perturbations,  respectively^^  In  the  latter  region,  the 
constant  term  is  irrelevant  and  V  may  be  approximated  by  —Vo  e-c<^  which 
may  be  studied  using  the  scaling  solution  discussed  in  Section  || 

For  an  arbitrary  scalar  potential  of  the  form  sketched,  i.e.  rising  with 
negative  curvature  towards  a  flat  plateau,  the  the  scalar  spectral  index  is 
given  approximately  by0li! 


ns  ~  1  —  4 


1  + 


V"V 

(V'Y 


(2) 


to  be  evaluated  when  the  modes  on  the  length  scales  of  interest  are  generated, 
Stage  (4)  as  described  in  Fig.  1.  For  the  exponential  form  here,  Eq.  © 
reduces  to 

ns  ~  1  4-  (3) 

Current  observational  limits  from  the  cosmic  microwave  background  and 
large  scale  structure  data  are  safely  satisfied  for  c  =  10,  which  we  shall 
adopt  as  our  canonical  value. 

The  fact  that  the  potential  minimum  is  negative  means  that  there  are 
no  strictly  static  solutions  for  0  except  anti-de  Sitter  space.  However,  as 
we  shall  show,  the  generic  behavior  -  indeed  an  attractor  -  is  a  dynamical 
“hovering”  solution  in  which  0  roams  back  and  forth  in  cyclic  fashion  between 
the  plateau  and  —  oo.  The  hovering  solution  is  highly  asymmetric  in  time. 
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Figure  1:  The  interbrane  potential  V((f>)  versus  (p,  whose  value  (—00  <  (p  < 
determines  the  distance  between  branes.  The  shaded  circle  represents 
the  maximum  positive  value  of  (f)  during  the  cycle.  The  various  stages  are: 
(1)  quintessence/potential  domination  and  cosmic  acceleration  (duration  > 
trillion  years);  (2)  (p  kinetic  energy  becomes  non- negligible,  decelerated  ex¬ 
pansion  begins  (duration  ~  1  billion  years);  (3)  H  =  0,  contraction  begins; 
(4)  density  fluctuations  on  observed  scales  created  ((toti?)1^2  ~  1  ms  before 
big  crunch);  (5)  (p  kinetic  energy  domination  begins  (fmin  ~  1CT30  s  before 
big  crunch);  (6)  bounce  and  reversal  from  big  crunch  to  big  bang;  (7)  end 
of  (p  kinetic  energy  domination,  potential  also  contributes  {tm/m  ~  1CT30  s 
after  big  bang);  (8)  radiation  dominated  epoch  begins  ~  1CT25  s  after  big 
bang);  (9)  matter  domination  epoch  begins  (~  1010  s  after  big  bang).  As  the 
potential  begins  to  dominate  and  the  Universe  returns  to  stage  (1),  the  field 
turns  around  and  rolls  back  towards  —00. 
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The  field  <f>  spends  trillions  of  years  or  more  on  the  plateau  and  mere  instants 
traveling  from  the  potential  well  to  —  oo  and  back.  Gravity  and  the  bounce 
provide  transfers  of  gravitational  to  kinetic  to  matter-radiation  density  that 
keep  the  Universe  forever  hovering  around  the  anti-de  Sitter  minimum  rather 
than  being  trapped  in  it. 

2.2  The  View  from  Effective  4d  Theory 

To  set  the  context  for  the  later  sections,  we  present  a  brief  tour  through  a 
single  cycle,  using  the  labels  in  Figure  1  as  the  mileposts.  Stage  1  represents 
the  present  epoch.  The  current  value  of  the  Hubble  parameter  is  Hq  =  (15 
billion  yr)_1.  We  are  presently  at  the  time  when  the  scalar  field  is  acting  as 
a  form  of  quintessence  in  which  its  potential  energy  has  begun  to  dominate 
over  matter  and  radiation.  Depending  on  the  specific  details  of  the  potential, 
the  field  (j)  may  have  already  reached  its  maximal  value  (grey  circle),  turned 
back,  and  begun  to  evolve  towards  negative  values.  If  not,  it  will  do  so  in 
the  near  future.  Because  the  slope  of  the  potential  is  very  small,  (j)  rolls  very 
slowly  in  the  negative  direction.  As  long  as  the  potential  energy  dominates, 
the  Universe  undergoes  exceedingly  slow  cosmic  acceleration  (compared  to 
inflationary  expansion),  roughly  doubling  in  size  every  Hq1  =  15  billion 
years.  If  the  acceleration  lasts  trillions  of  years  or  more  (an  easy  constraint 
to  satisfy),  the  entropy  and  black  hole  densities  become  negligibly  small  and 
the  Universe  is  nearly  vacuous.  The  Einstein  equations  become: 


H 2  = 

87 tG  1 
3  \ 

'\i?  +  rw) 

(4) 

a 

8nG 

02  -  VOW) 

(5) 

a 

3 

where  H  is  the  Hubble  parameter  and  G  is  Newton’s  constant.  We  will  gen¬ 
erally  choose  87rG  =  1  except  where  otherwise  noted.  Accelerated  expansion 
stops  as  V((f>)  approaches  zero  and  the  scalar  field  kinetic  energy  becomes 
comparable  to  the  potential  energy,  Stage  (2).  The  Universe  continues  to 
expand  and  the  kinetic  energy  of  scalar  field  continues  to  red  shift  as  the 
potential  drops  below  zero.  A  nearly  scale  invariant  spectrum  of  fluctuations 
on  large  length  scales  (beyond  our  current  Hubble  horizon)  begins  to  develop 
as  the  field  rolls  down  the  exponentially  decreasing  part  of  the  potential. 
The  evolution  and  perturbation  equations  are  the  same  as  in  the  ekpyrotic 
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model .IZlliZl  Solving  these  equations,  one  finds  that  the  decelerated  expansion 
continues  for  a  time  Hq1/c:  which  is  about  one  billion  years  (c  is  the  pa¬ 
rameter  in  U(0),  Eq.  ([]])).  At  Stage  (3),  the  potential  becomes  sufficiently 
negative  that  the  total  scalar  field  energy  density  hits  zero.  According  to 
Eq.  ((H),  H  =  0  and  the  Universe  is  momentarily  static.  From  Eq.  (^),  a  <  0, 
so  that  a  begins  to  contract.  The  Universe  continues  to  satisfy  the  ekpyrotic 
conditions  for  creating  density  perturbations.  Stage  (4),  about  one  second 
before  the  big  crunch,  is  the  regime  where  fluctuations  on  the  current  Hubble 
horizon  scale  are  generated.  As  the  field  continues  to  roll  towards  —  oo,  the 
scale  factor  a  contracts  and  the  kinetic  energy  of  the  scalar  field  grows.  That 
is,  gravitational  energy  is  converted  to  scalar  field  (brane)  kinetic  energy  dur¬ 
ing  this  part  of  the  cycle.  Hence,  the  field  races  past  the  minimum  of  the 
potential  at  Stage  (5)  and  off  to  — oo,  with  kinetic  energy  becoming  increas¬ 
ingly  dominant  as  the  bounce  approaches.  The  scalar  field  kinetic  energy 
diverges  as  a  tends  to  zero.  At  the  bounce,  Stage  (6),  matter  and  radiation 
are  generated,  the  scalar  field  gets  a  kick  and  increases  speed  as  it  reverses 
direction,  and  the  Universe  is  expanding.  Through  Stage  (7),  the  scalar  ki¬ 
netic  energy  density  (oc  1/a6)  dominates  over  the  radiation  (oc  1/a4)  and  the 
motion  is  almost  exactly  the  time-reverse  of  the  contraction  phase  between 
Stage  (5)  and  the  big  crunch.  As  the  field  rolls  uphill,  however,  the  small 
kick  given  the  scalar  field  and,  subsequently,  the  radiation  become  impor¬ 
tant,  breaking  the  time-reversal  symmetry.  The  Universe  becomes  radiation 
dominated  at  Stage  (8),  at  say  10-25  s  after  the  big  bang.  The  motion  of  0  is 
rapidly  damped  so  that  it  converges  towards  its  maximal  value  and  then  very 
slowly  creeps  downhill.  The  damping  continues  during  the  matter  dominated 
phase,  which  begins  thousands  of  years  later.  The  Universe  undergoes  the 
standard  big  bang  evolution  for  the  next  15  billion  years,  growing  structure 
from  the  perturbations  created  when  the  scalar  field  was  rolling  downhill  at 
Stage  (4).  Then,  the  scalar  field  potential  energy  begins  to  dominate  and  cos¬ 
mic  acceleration  begins.  Eventually,  the  scalar  field  rolls  back  across  0  =  0. 
The  energy  density  falls  to  zero  and  cosmic  contraction  begins.  The  scalar 
field  rolls  down  the  hill,  density  perturbations  are  generated  and  0  runs  off 
to  — oo  for  the  next  bounce.  The  evolution  in  terms  of  conventional  variables 
is  summarized  in  Figure  [2.2|. 
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log  t 

Figure  2:  Schematic  plot  of  the  scale  factor  a(t),  the  modulus  </>(£),  and 
H5  =  | d (exp ( y^3/2 (j))/dt  for  one  cycle,  where  t  is  Einstein  frame  proper 

time.  The  scale  factor  starts  out  zero  but  expands  as  1 3,  and  the  scalar 

held  grows  logarithmically  with  f,  in  the  scalar  kinetic  energy  dominated 

early  regime.  Then,  when  radiation  begins  to  dominate  we  have  a  oc  and 

the  scalar  held  motion  is  strongly  damped.  This  is  followed  by  the  matter 
2 

era,  where  a  (x  1 3 ,  and  a  potential  dominated  phase  in  which  a(t)  increases 
exponentially,  before  a  final  collapse  on  a  timescale  Hq1,  to  a  =  0  once  more. 
H5  is  proportional  to  the  proper  (five  dimensional)  speed  of  contraction  of 
the  fifth  dimension.  To  obtain  a  cyclic  solution,  the  magnitude  of  H5  at  the 
start  of  the  big  bang,  H5(out),  must  be  slightly  larger  than  the  value  at  the 
end  of  the  big  crunch,  H5(in).  This  is  the  case  if  more  radiation  is  generated 
on  the  negative  tension  brane  (see  Appendix). 
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2.3  The  View  from  the  Visible  Brane 


Thus  far,  we  have  described  the  evolution  in  terms  of  the  usual  Einstein 
frame  variables,  a  and  0.  However,  as  emphasized  in  the  next  Section,  these 
variables  are  singular  at  the  transition  from  big  bang  to  big  crunch,  and  they 
do  not  present  an  accurate  picture  of  what  an  observer  composed  of  matter 
confined  to  the  brane  would  actually  see.  As  a  approaches  zero,  the  density  of 
matter  and  radiation  scale  as  l/(a0(0))3  and  l/(a0(0))4,  respectively,  where 
0(0)  is  a  function  of  0  which  scales  as  1/a  as  a  tends  to  zero.  Therefore  the 
densities  of  matter  and  radiation  on  the  branes  are  actually  finite  at  a  =  0. 

This  scaling  of  the  density  with  a0(0)  rather  than  a  can  be  understood 
rather  simply.  First,  the  spatial  volume  element  on  the  branes  is  that  induced 
from  five  dimensions.  When  the  brane  separation  is  small,  one  can  use  the 
usual  formula  for  Kaluza-Klein  theory, 

dsl  =  e“V^ds2  +  e2yf^<t>dy2,  (6) 

where  ds\  is  the  four  dimensional  line  element,  y  is  the  fifth  spatial  coordinate 
which  runs  from  zero  to  L ,  and  L  is  a  parameter  with  the  dimensions  of 
length.  If  we  write  the  four  dimensional  line  element  in  conformal  time 
coordinates,  as  ds\  =  a2(— d,T2-{-dx2),  then  since  from  the  Friedmann  equation 
we  have  ( a' /a )2  =  |(0')2,  we  see  that  a  is  proportional  to  e in  the  big 

crunch.  Hence  a  three  dimensional  comoving  volume  element  d3xa3e~ 
remains  finite  as  a  tends  to  zero.  Thus  the  density  of  massive  particles  tends 
to  a  constant.  What  about  the  density  of  radiation?  First,  recall  the  usual 
argument  that  the  energy  of  a  photon  diverges  at  a  =  0.  Consider  a  set  of 
comoving  massive  particles  in  a  spacetime  with  metric  a 2r)lxv  where  r)lw  is  the 
Minkowski  metric.  The  four  velocities  of  the  particles  obey  u,luu  gliv  =  —  1. 
Hence,  if  they  are  comoving  (u  —  0),  then  we  must  have  u°  =  a-1.  Now 
a  photon  moving  in  such  a  spacetime  has  a  constant  four-momentum,  p M  = 
E(l,n),  with  n 2  =  1.  The  energy  of  the  photon,  as  seen  by  the  comoving 
particles,  is  —uflpugfll/  =  E/a,  which  diverges  a  tends  to  zero.  However, 
in  the  present  context,  the  metric  to  which  the  comoving  particles  couple 
is  e^V^a2//^.  Therefore,  we  have  u°  =  a-1eV^  and  the  energy  of  the 
detected  photons  is  finite  as  a  tends  to  zero.  In  other  words,  the  effect  of  the 
scalar  field  approaching  —  oo  is  precisely  such  as  to  cancel  the  gravitational 
blueshift. 
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The  second  crucial  use  of  the  higher  dimensional  metric  is  in  piecing  to¬ 
gether  the  global  view  of  the  spacetime.  If  one  only  had  the  Einstein  frame 
scale  factor  a,  it  would  not  be  clear  how  to  match  to  the  next  cycle,  since 
a  =  0  at  the  bounce.  But  the  scale  factor  on  a  brane,  a/3(4>),  is  nonzero  at 
each  bounce  and  may  be  so  matched.  In  fact,  in  the  examples  studied  in  this 
paper,  the  scale  factors  o0  and  a\  (which  are  the  brane  scale  factors  in  the 
simplest  models)  both  undergo  a  net  exponential  expansion  within  a  cycle, 
and  decrease  for  very  brief  periods  -  either  just  before  the  brane  collision  (for 
«o)  or  just  after  it  (for  op).  An  observer  on  either  brane  would  view  the  cos¬ 
mology  as  one  of  almost  uninterrupted  expansion,  with  successive  episodes  of 
radiation-,  matter-,  and  quintessence-domination  ending  in  a  sudden  release 
of  matter  and  radiation. 

Both  matter  and  radiation  are  suddenly  created  by  the  impact  of  the 
other  brane.  The  forewarning  of  this  catastrophic  event  would  be  that  as 
/3(0)  started  to  rapidly  change,  one  would  see  stronger  and  stronger  violations 
of  the  equivalence  principle  (a  ‘fifth  force’),  and  the  masses  and  couplings  of 
all  particles  would  change.  In  the  case  of  M  theory,  the  running  of  the  string 
coupling  to  zero  would  presumably  destroy  all  bound  states  such  as  nucleons 
and  send  all  particle  masses  to  zero. 

3  Through  the  Bounce 

To  have  repeating  cycles,  the  Universe  must  be  able  to  pass  smoothly  from 
a  big  crunch  to  a  big  bang.  Conventionally,  the  curvature  and  density  sin¬ 
gularity  when  the  scale  factor  a  approaches  zero  has  been  regarded  as  an 
impassable  obstacle  to  the  understanding  of  what  came  ‘before’  the  big  bang. 
However,  the  brane  world  setup  sheds  new  light  on  this  problem.  The  key 
feature  is  that  the  apparent  singularity  in  the  effective  four-dimensional  de¬ 
scription  corresponds  to  a  higher  dimensional  setup  in  which  the  four  dimen¬ 
sional  metric  is  completely  non-singular.  When  the  extra  dimension  (or  outer 
brane  separation)  shrinks  to  zero,  there  is  no  associated  curvature  singularity, 
and  the  density  of  matter  on  the  branes  remains  finite.  The  most  conserva¬ 
tive  assumption,  based  on  the  higher  dimensional  picture,  is  that  the  branes 
bounce  from  (or,  equivalently,  pass  through)  each  other  and  time  continues 
smoothly,  with  some  conversion  of  brane  kinetic  energy  to  entropy.  The  sep¬ 
aration  of  the  two  branes  after  the  bounce  corresponds  to  re-expansion  in 
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the  four- dimensional  effective  theory. 

How  can  this  be  reconciled  with  the  singular  four-dimensional  descrip¬ 
tion?  The  point  explained  in  Ref.  |8]  is  that  the  usual  four-dimensional  vari¬ 
ables,  the  scale  factor  a  and  the  scalar  field  (ft ,  are  a  singular  choice  at  a  =  0. 
Each  is  poorly  behaved  as  the  branes  collide,  but  in  the  brane  picture  physi¬ 
cal  quantities  depend  on  combinations  of  the  two  variables  that  remain  well- 
behaved.  These  nonsingular  variables  may  be  treated  as  fundamental,  and 
matching  rules  derived  to  parameterize  the  physics  of  inelastic  brane  colli¬ 
sions.  If  the  system  can,  as  conjectured  in  Ref.  |],  be  properly  embedded 
within  string  theory,  the  matching  conditions  will  be  derivable  from  funda¬ 
mental  physics. 


3.1  Non-singular  variables 

The  action  for  a  scalar  held  coupled  to  gravity  and  a  set  of  fluids  pi  in  a 
homogeneous,  hat  Universe,  with  line  element  ds2  =  a2(r)(—N2dT2  +  dx2)  is 


5 


d3xdr 


-N((a(3)4XiPi  +  a4V((ft))) 


(7) 


We  use  r  to  represent  conformal  time  and  primes  to  represent  derivatives 
with  respect  to  r.  N  is  the  lapse  function.  The  background  solution  for  the 
scalar  held  is  denoted  0(r),  and  V(<p)  is  the  scalar  potential. 

The  only  unusual  term  in  (0)  is  the  coupling  of  the  fluids  pt,  which  we 
treat  as  perfect  fluids  coupled  only  through  gravity.  The  action  for  a  perfect 
fluid  coupled  to  gravity  is  just  —  /  d4x^/—gp,  where  the  density  p  is  regarded 
as  a  function  of  the  coordinates  of  the  fluid  particles  and  the  spacetime  met- 
ric.0  For  a  homogeneous  isotropic  fluid,  the  equation  of  state  P(p)  dehnes 
the  functional  dependence  of  p  on  the  scale  factor  a,  via  energy-momentum 
conservation,  dlnp/dlna  =  —3(1  +  w),  with  w  =  P/p.  For  example,  for 
radiation,  p  oc  a-4  and  for  matter  p  oc  cT3. 

We  assume  these  fluids  live  on  one  of  the  branes,  so  that  rather  than 
coupling  to  the  Einstein-frame  scale  factor  a,  the  particles  they  are  composed 
of  couple  to  a  conformally  related  scale  factor  a/3(0),  being  the  scale  factor 
on  the  appropriate  brane.  For  simplicity  we  have  only  written  the  action  for 
fluids  on  one  of  the  branes,  the  action  for  fluids  on  the  other  brane  being  a 
xerox  copy  but  with  the  appropriate  /3(cft). 
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The  function  0(0)  may  generally  be  different  for  the  two  branes,  and  for 
different  brane  world  setups.  But  as  mentioned  above  there  is  an  important 
universality  at  small  separations  corresponding  to  large  negative  0.  In  this 
limit,  which  is  relevant  to  the  bounce,  the  bulk  warp  factor  becomes  irrelevant 
and  one  obtains  0  ~  e-9Vv6,  the  standard  Kaluza-Klein  result.  This  behavior 
ensures  that  a/3  is  finite  at  collision  and  so  the  matter  and  radiation  densities 
are,  as  well. 

The  equations  of  motion  for  gravity,  the  matter  and  scalar  field  0  are 
straightforwardly  derived  by  varying  (0)  with  respect  to  a,  N  and  0,  after 
which  N  may  be  set  equal  to  unity.  Expressed  in  terms  of  proper  time  t,  The 
Einstein  equations  are 


H2  = 


87 tG 


1  ; 
2 


0“  +  V  +  (34Pr  +  (34pM^j  , 


(8) 


; =  -  ■ if  02  - v + p‘Pr + \i3‘Pm)  •  (9) 

where  a  dot  is  a  proper  time  derivative.  As  an  example,  we  consider  the 
case  where  there  is  radiation  (pr)  and  matter  (pm)  on  the  visible  brane  only, 
which  could  in  principle  be  either  the  positive  or  negative  tension  brane. 
Then  the  above  equations  are  supplemented  by  the  dynamical  equation  for 
the  evolution  of  0, 

0  +  3H<j>  =  -  (3^pM  (10) 

and  the  continuity  equation, 


^dp  dp  [3  dp 

a  da,  da  +  f3'  d<f) 


-3  (p  +  p) 


(11) 


where  a  =  af3{4>)  and  p  is  the  pressure  of  the  fluid  component  with  energy 
density  p.  Note  that  only  the  matter  density  contributes  to  the  0-equation, 
because,  if  pr  oc  1  / (a/3)4,  the  radiation  term  is  actually  just  a  constant 
times  N  in  the  action,  contributing  to  the  Friedmann  constraint  but  not  the 
dynamical  equations  of  motion. 

If  0(0)  is  sufficiently  flat  near  the  current  value  of  0,  these  couplings 
have  modest  effects  in  the  late  Universe,  and  the  successes  of  the  standard 
cosmology  are  recovered.  For  example  the  total  variation  in  0  since  nucle¬ 
osynthesis  is  very  modest.  I11  Planck  units,  this  is  of  order  (tr/tjv)^  where 
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tr  is  the  time  at  the  beginning  of  the  radiation  dominated  epoch  and  nu¬ 
cleosynthesis  begins  at  £/v  ~  1  sec.  It  is  utterly  negligible  for  values  of  tr 
earlier  than  the  electroweak  era.  However,  the  coupling  of  matter  to  0  pro¬ 
duces  other  potentially  measurable  effects  including  a  ‘fifth  force’  causing 
violations  of  the  equivalence  principle.  Current  constraints  can  be  satisfied 
if  MPl(ln/3)t4>  <  10“3.&E3 

As  the  Universe  contracts  towards  the  big  crunch,  a  — »  0,  the  scalar  held 
runs  to  — oo  and  the  scalar  potential  becomes  negligible.  The  Universe  be¬ 
comes  dominated  by  the  scalar  held  kinetic  energy  density  since  it  scales  as 
a-6  whereas  matter  and  radiation  densities  scale  as  a-3  and  a-4  respectively 
(ignoring  the  f3  factor).  As  scalar  kinetic  domination  occurs,  the  scale  factor 
a  begins  to  scale  as  (—  f)s,  and  the  background  scalar  held  diverges  logarith¬ 
mically  in  time.  The  energy  density  and  Ricci  scalar  diverge  as  (— 1)~2,  so 
that  t  =  0  is  a  ‘big  crunch’  singularity. 

As  explained  in  Ref.  in  the  simplest  treatment  of  brane  world  models 
there  is  only  one  scalar  held  modulus,  the  ‘radion,’  which  runs  oh  to  minus 
infinity  as  the  scale  factor  a  approaches  zero.  The  singular  variables,  a  and 
0,  can  be  replaced  by  the  non-singular  variables: 

a0  =  2acosh((0  —  0oo)/V6)  ai  =  — 2  asinh((0  —  0oo)/a/6).  (12) 

The  kinetic  terms  in  the  action  dehne  the  metric  on  moduli  space.  In  terms 
of  the  old  variables  one  has  the  line  element  —3 da2  +  |a2d02,  and  a  =  0  is 
clearly  a  singular  point  in  these  coordinates.  However,  in  the  new  coordinates 
in  Eq.  (|12|),  the  line  element  is  |(— da^  +  daf),  which  is  perfectly  regular  even 

when  the  Einstein  frame  scale  factor  a  =  \\j  Oq  —  a2  vanishes,  on  the  ‘light- 
cone’  ao  =  ai.  For  branes  in  AdS,  a0  and  a\  are  the  scale  factors  on  the 
positive  and  negative  tension  branes3  so  that  (5  =  2cosh((0  —  0oo)/v/6  or 
— 2sinh((0  —  0oo)/v/6)  respectively  for  matter  coupling  to  these  branes. 

Notice  that  the  constant  held  shift  0^  is  arbitrary.  Its  effect  is  a  Lorentz 
boost  on  the  (a0,ai)  moduli  space.  In  the  Kaluza-Klein  picture  (|6]),  a  con¬ 
stant  shift  in  0  can  be  removed  by  rescaling  four  dimensional  spacetime 
coordinates  and  re-defining  the  length  scale  L  of  the  extra  dimensions.  In 
the  absence  of  matter  which  couples  to  0,  or  of  a  potential  U(0),  this  shift  is 
unobservable,  a  reflection  of  the  global  symmetry  0  — >  0+  constant  of  the  4d 
effective  theory.  However,  this  symmetry  is  broken  by  U(0),  and  by  matter 
couplings.  In  fact,  the  scale  factor  a\  must  be  positive  in  order  for  it  to  be 
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interpretable  as  a  ‘brane  scale  factor’,  and  this  requires  that  0  <  0oc- 

We  shall  find  it  convenient  to  choose  0  =  0  to  be  the  zero  of  the  potential 
V (0) ,  and  then  to  choose  0^  so  that  <q  never  vanishes  for  the  solutions  we 
are  interested  in.  (In  fact,  since  cq  has  a  positive  kinetic  terrn,§  a  suitable 
coupling  to  moduli  fields  will  always  guarantee  that  cq  ‘bounces’  away  from 
zeroJil  In  this  paper,  for  simplicity  we  ignore  this  complication  by  picking 
0oo  large  enough  that  no  such  ‘bounce’  is  necessary.) 

Both  ao  and  cq  are  ‘scale  factors’  since  they  transform  like  a  under  rescal¬ 
ing  space-time  coordinates.  However,  unlike  a  they  tend  to  finite  constants 
as  a  tends  to  zero,  implying  an  alternative  metrical  description  which  is  not 
singular  at  the  ‘big  crunch’.  In  the  brane  world  models  considered  in  Ref.  [F|. 
cio  and  d\  actually  represent  the  scale  factors  of  the  positive  and  negative  ten¬ 
sion  branes  respectively.  Since  there  are  no  low  energy  configurations  with 
a0  <  di,  the  ‘light  cone’  a0  =  cq  is  actually  a  boundary  of  moduli  space  and 
one  requires  a  matching  rule  to  determine  what  the  trajectory  of  the  system 
does  at  that  point.  A  natural  matching  rule  is  to  suppose  that  at  low  energies 
and  in  the  absence  of  potentials  or  matter,  the  branes  simply  pass  through 
one  another  (or,  equivalently,  bounce)  with  the  intervening  bulk  briefly  dis¬ 
appearing  and  then  reappearing  after  collision.  This  rule  was  detailed  in 
Ref.  where  simple  models  satisfying  the  string  theory  background  equa¬ 
tions  to  all  orders  in  a'  were  given.  In  the  Appendix  we  discuss  the  collision 
between  boundary  branes  in  terms  of  energy  and  momentum  conservation, 
and  the  Israel  matching  conditions. 

Let  us  now  comment  on  the  character  of  the  trajectory  in  the  (do,  cq)- 
plane.  The  Friedmann  constraint  reads 

«a  ~a'i  =  l  ((^)4P+^(a  o  -  o))  •  (13) 

If  the  energy  density  on  the  right  hand  side  is  positive,  the  trajectory  is 
time-like.  If  the  right  hand  side  is  zero  (for  example  if  the  potential  vanishes 
as  0o  — * ►  —  oo  and  if  there  is  no  matter  or  radiation),  then  the  trajectory  is 
light-like.  If  the  right  hand  side  is  negative,  the  trajectory  is  space-like. 

The  trajectory  for  the  cyclic  solution  in  the  a0  —  eq  plane  is  shown  in 
Figure  |Q|.  The  insert  shows  a  blow-up  of  the  behavior  at  the  bounce  in  which 
the  trajectory  is  light-like  at  contraction  to  the  big  crunch  (the  Universe  is 
empty)  and  time-like  on  expansion  from  the  big  bang  (radiation  is  produced 
at  the  bounce).  In  these  coordinates,  the  scale  factor  increases  exponentially 
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over  each  cycle,  but  the  next  cycle  is  simply  a  rescaled  version  of  the  cycle 
before.  A  local  observer  measures  physical  quantities  such  as  the  Hubble 
constant  or  the  deceleration  parameter,  which  entail  ratios  of  the  scale  factor 
and  its  derivatives  in  which  the  normalization  of  the  scale  factor  cancels  out. 
Hence,  to  local  observers,  each  cycle  appears  to  be  identical  to  the  one  before. 

3.2  From  Big  Crunch  to  Big  Bang 

In  this  section  we  solve  the  equations  of  motion  immediately  before  and  after 
the  bounce,  and  discuss  how  the  incoming  and  outgoing  states  are  connected. 
The  nonsingular  ‘brane’  scale  factors  «o  and  cq  provide  the  natural  setting 
for  this  discussion,  since  neither  vanishes  at  the  bounce.  As  emphasized 
above,  the  Einstein  frame  scale  factor  a,  and  the  scalar  field  0  are  singular 
coordinates  on  field  space  at  the  bounce.  Nevertheless,  since  onr  intuition  is 
much  better  in  the  Einstein  frame,  we  shall  also  give  formulae  for  a  and  0 
near  the  bounce.  In  subsequent  sections,  we  shall  frame  the  discussion  almost 
entirely  in  terms  of  Einstein  frame  variables,  for  the  most  part  using  the 
nonsingular  variables  «o  and  ci\  solely  as  a  ‘bridge’  connecting  the  incoming 
big  crunch  to  the  outgoing  big  bang. 

Before  the  bounce  there  is  little  radiation  present  since  it  has  been  expo¬ 
nentially  diluted  in  the  preceding  quintessence-dominated  accelerating  phase. 
Furthermore,  the  potential  V(<p)  becomes  negligible  as  0  runs  off  to  minus 
infinity.  The  Friedmann  constraint  reads  (a'/a)2  =  |0/2,  and  the  scalar  field 
equation,  (a20/)/  =  0,  where  primes  denote  conformal  time  derivatives.  The 
general  solution  is 

0  =  ^^ln  (AH5(in)r) ,  a  =  =  A\J AH5(in)r, 

ao  =  A  (^X  +  X~1AH5(in)TSj,  a\  —  A  (A  —  A-1  AH^iin)^  ,  (14) 

where  A  =  .  We  choose  r  =  0  to  be  the  time  when  a  vanishes  so 

that  r  <  0  before  collision.  A  is  an  integration  constant  which  could  be  set 
to  unity  by  rescaling  space-time  coordinates  but  it  is  convenient  not  to  do 
so.  The  Hubble  constants  as  defined  in  terms  of  the  brane  scale  factors  are 
«q/ Oq  and  a) / a2  which  at  r  =  0  take  the  values  +X~3H5(in)  and  —X~3H5(in) 
respectively. 
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Figure  3:  Schematic  plot  of  the  a0-Gq  plane  showing  a  sequence  of  cycles 
of  expansion  and  contraction  (indicated  by  tick  marks).  The  dashed  line 
represents  the  “light-cone”  a0  =  cq  corresponding  to  a  bounce  (a  =  0).  Each 
cycle  includes  a  moduli  kinetic  energy,  radiation,  matter  and  quintessence 
dominated  phase  and  lasts  an  exponentially  large  number  of  e-folds.  The 
insert  shows  the  trajectory  near  the  big  crunch  and  bounce.  The  potential 
energy  V{<p)  assumed  takes  the  form  shown  in  Fig.  1. 
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Re-expressing  the  scalar  field  as  a  function  of  proper  time  t  —  f  adr,  we 
obtain 


</>  =  yfin  . 


(15) 


The  integration  constant  H^{in)  <  0  has  a  natural  physical  interpretation  as 
a  measure  of  the  contraction  rate  of  the  extra  dimension  (See  equation  (^|)): 


H,  = 


dL5  _  d  (e v  3 


?0e\/T* 


(16) 


Ldt5  dt5 

where  L5  =  Le VlN’  is  the  proper  length  of  the  extra  dimension,  L  is  a 
parameter  with  dimensions  of  length,  and  1 5  is  the  proper  time  in  the  five¬ 
dimensional  metric, 

dt5  =  ae~^dr  =  e~^dt,  (17) 

with  t  being  FRW  proper  time.  Notice  that  a  shift  (j) ^  can  always  be  com¬ 
pensated  for  by  a  rescaling  of  L.  As  the  extra  dimension  shrinks  to  zero, 
tends  to  a  constant,  H5(in). 

If  the  extra  dimension  shrinks  adiabatically  and  backreaction  from  parti¬ 
cle  production  can  be  ignored,  then  the  matching  rule  conjectured  in  Ref.  |8] 
states  that  H5  after  the  bounce  should  be  given  by  H5(out )  = 

However,  if  radiation  is  produced,  H5(out)  takes  a  different  value.  If  one 
is  given  the  densities  of  radiation  produced  on  both  branes,  then  H5(out ) 
may  be  inferred  from  energy  and  momentum  conservation,  and  the  Israel 
matching  conditions,  as  we  show  in  the  Appendix. 

Immediately  after  the  bounce,  scalar  kinetic  energy  dominates  and  Hr, 
remains  nearly  constant,  as  shown  in  Figure  |2.2|.  The  kinetic  energy  of  the 
scalar  field  scales  as  a-6  and  radiation  scales  as  a-4,  so  the  former  dominates 
at  small  a.  It  is  convenient  to  re-scale  a  so  that  it  is  unity  at  scalar  kinetic 
energy-radiation  equality,  tr,  and  denote  the  corresponding  Hubble  constant 
Hr.  The  Friedmann  constraint  in  Eq.  (13)  then  reads 


(«/)2=2^2(1  +  °-2)’ 


and  the  solution  is 


2  a  rH^  (out) H? 
(Hrr  +  2§) 


a  =  J  —H?t2  +  V2Ht.t. 


(18) 


(19) 
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The  brane  scale  factors  are 


a0  =  a(A"1e^/V®  +  \e~^)  =  A  (a(1  +  ^f)  +  A H? H5§ (out)r')  , 
ai  =  a(-A~1e(/,/^  +  Ae“*/V®) 

=  A  ^A(l  +  )  -  2*A~1H}Hl(out)TSj  .  (20) 

1  1 

Here  the  constant  A  —  2&  (Hr/ H^(out))3  has  been  defined  so  that  we  match 
ao  and  (i\  to  the  incoming  solution  given  in  (|T|).  As  for  the  incoming  solution, 

we  can  compute  the  Hubble  constants  on  the  two  branes  after  collision.  They 

5  2  1 

are  ±\~3H5(out)  +  2~3  \~lHr  on  the  positive  and  negative  tension  branes 

respectively. 

For  Hr  <  22\~3H5,  the  case  of  relatively  little  radiation  production, 
immediately  after  collision  a 0  is  expanding  but  cq  is  contracting.  Whereas 
for  Hr  >  2a A  3H5,  both  brane  scale  factors  expand  after  collision.  We 
shall  concentrate  on  the  former  case  in  this  paper,  in  which  we  are  near  the 
adiabatic  limit.  If  no  scalar  potential  V((p)  were  present,  the  scalar  field 
would  continue  to  obey  the  solution  (®.  converging  to 



(21) 


This  value  is  actually  larger  than  0^  for  Hr  <  H5 A-325,  the  case  of  weak 
production  of  radiation.  However,  the  presence  of  the  potential  H(0)  alters 
the  expression  (2l|)  for  the  final  resting  value  of  the  scalar  field.  As  0  crosses 
the  potential  well  traveling  in  the  positive  direction,  H5  is  reduced  to  a 
renormalized  value  H5(out )  <  H5(out),  so  that  the  final  resting  value  of  the 
scalar  field  can  be  smaller  than  0^.  If  this  is  the  case,  then  cp  never  crosses 
zero,  instead  reversing  to  expansion  shortly  after  radiation  dominance.  (In 
the  calculations  of  Ref  O,  where  we  assumed  the  potential  vanished  after 


collision,  this  effect  did  not  occur.  Instead,  we  invoked  a  coupling  of  cp  to  a 
modulus  field  which  caused  it  to  bounce  off  cp  =  0.). 

If  radiation  dominance  occurs  well  after  0  has  crossed  the  potential  well, 
Eq.  (|2l|)  provides  a  reasonable  estimate  for  the  final  resting  value,  if  we  use 
the  corrected  value  H^out).  The  dependence  of  (f2T|)  is  simply  understood: 
while  the  Universe  is  kinetic  energy  dominated,  a  grows  as  fa  and  0  increases 
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logarithmically  with  time.  However,  when  the  Universe  becomes  radiation- 
dominated  and  a  oc  £2,  Hubble  damping  increases  and  0  converges  to  the 
finite  limit  above. 


4  Across  the  Well 


Using  the  potential  described  in  Section  ^T]  and,  specifically,  the  example 
in  Eq.  fll]),  this  Section  considers  the  motion  of  0  back  and  forth  across  the 
potential  well.  We  will  show  that  evolution  converges  to  a  stable  attractor 
solution.  Our  main  purpose,  though,  is  to  explore  the  asymmetry  in  the 
behavior  before  and  after  the  bounce  that  is  an  essential  component  of  the 
cyclic  solution. 

Over  most  of  this  region,  V  may  be  accurately  approximated  by  —  Vn  ■ 
For  this  pure  exponential  potential,  there  is  a  simple  scaling  solution&3 


a(t)  =  \t\p, 


V  =  -V0e-^  = 


P( 1  -  3 p) 
t2  ' 


P  = 


(22) 


which  is  an  expanding  or  contracting  Universe  solution  according  to  whether 
£  is  positive  or  negative.  (We  choose  £  =  0  to  be  the  bounce.)  From  the 
expression  for  V,  we  see  that  0  varies  logarithmically  with  \t\. 

At  the  end  of  the  expanding  phase  of  the  cyclic  scenario,  there  is  a  period 
of  accelerated  expansion  which  makes  the  Universe  empty,  homogeneous  and 
flat,  followed  by  0  rolling  down  the  potential  U(0)  into  the  well.  After  0  has 
rolled  sufficiently  and  the  scale  factor  has  begun  to  contract  (past  Stage  (3)  in 
Figure  |0|),  the  Universe  accurately  follows  the  above  scaling  solution  down 
the  well  until  0  encounters  the  potential  minimum  (Stage  (5)  in  Figure  |2TT|) . 

Let  us  consider  the  behavior  of  0  under  small  shifts  in  the  contracting 
phase.  In  the  background  scalar  field  equation  and  the  Friedmann  equation, 
we  set  0  =  4>b  +  <50  and  H  =  HB  +  5H,  where  0^  and  HB  are  the  background 
quantities  given  from  (|22|).  To  linear  order  in  50,  one  obtains 


50  +  ^^50 


1  -  3  p 


50  =  0, 


(23) 


with  two  linearly  independent  solutions,  50  ~  t~l  and  f1-3p,  where  p  • Cl.  In 
the  contracting  phase,  the  former  solution  grows  as  t  tends  to  zero.  However, 
this  solution  is  simply  an  infinitesimal  shift  in  the  time  to  the  big  crunch: 
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<50  oc  0.  Such  a  shift  provides  a  solution  to  the  Einstein-scalar  equations 
because  they  are  time  translation  invariant,  but  it  is  physically  irrelevant 
since  it  can  be  removed  by  a  redefinition  of  time.  The  second  solution  is  a 
physical  perturbation  mode  and  it  decays  as  t  tends  to  zero.  Hence,  we  find 
that  the  background  solution  is  an  attractor  in  the  contracting  phase. 

We  next  consider  the  incoming  and  outgoing  collision  velocity,  which  we 
have  parameterized  as  H5(in )  and  H$(out)  in  the  previous  Section.  Within 
the  scaling  solution  (p2|),  we  can  calculate  the  value  of  incoming  velocity  by 
treating  the  prefactor  of  the  potential  F(0)  in  Eq.  (!l|)  as  a  Heaviside  function 
which  is  unity  for  0  >  0min  and  zero  for  0  <  0min,  where  <pmin  is  the  value  of 
0  at  the  minimum  of  the  potential.  We  compute  the  velocity  of  the  field  as 
it  approaches  (f)min  and  use  energy  conservation  at  the  jump  in  V  to  infer  the 
velocity  after  0m,n  is  crossed.  In  the  scaling  solution,  the  total  energy  as  0 
approaches  0min  from  the  right  is  |02  +  V  =  3 p2/t2,  and  this  must  equal  the 
total  energy  |02  evaluated  for  0  just  to  the  left  of  0min .  Hence,  we  find  that 
0  =  \/6 p/t  =  yj QpVmin/ (1  —  3 p)  at  the  minimum  and,  according  to  Eq.  (jT6|), 


H5(in )  « 


V8\Vmin\h^ 

c  y/1  —  6c-2 


(24) 


At  the  bounce,  this  solution  is  matched  to  an  expanding  solution  with 


H5(out )  =  -(1  +  x)H5(in)  >  0, 


(25) 


where  y  is  a  small  parameter  which  arises  because  of  the  inelasticity  of  the 
collision. 

In  order  to  obtain  cyclic  behavior,  we  shall  need  y  to  be  positive  or, 
equivalently,  the  outgoing  velocity  to  exceed  the  incoming  velocity.  There 
are  at  least  two  effects  that  can  cause  y  to  be  positive.  First,  as  we  discuss 
in  the  Appendix,  y  is  generically  positive  if  more  radiation  is  generated  on 
the  negative  tension  brane  than  on  the  positive  tension  brane  at  collision. 
Secondly,  y  can  get  a  positive  contribution  from  the  coupling  of  0(0)  to  the 
matter  created  on  the  branes  by  the  collision;  see  Eq.  (0).  Both  effects  are 
equally  good  for  our  purposes.  For  the  present  discussion,  we  shall  simply 
assume  a  small  positive  y  is  given,  and  follow  the  evolution  forwards  in  time. 

Since  y  is  small,  the  outgoing  solution  is  very  nearly  the  time  reverse 
of  the  incoming  solution  as  0  starts  back  across  the  potential  well  after  the 


25 


bounce:  the  scaling  solution  given  in  (|22|),  but  with  t  positive.  As  time 
proceeds  however,  the  contribution  of  x  becomes  increasingly  significant.  In 
the  time-reversed  scaling  solution,  H5  tends  to  zero.  For  x  >  0  H5  remains 
positive  and  0  overshoots  the  potential  well.  Vo  is  exponentially  smaller  than 
the  kinetic  energy  density  at  the  bounce,  so  even  a  tiny  fraction  x  suffices  to 
reach  the  plateau  after  crossing  the  potential  well. 

We  can  analyze  this  overshoot  by  treating  x  as  a  perturbation  and  using 
the  solution  in  Eq.  (|3|)  discussed  above,  5(ft  ~  t~1  and  f1_3p.  The  latter  is 
a  decaying  mode  in  the  contracting  phase  before  the  bounce  but  it  grows  in 
the  expanding  phase.  One  can  straightforwardly  compute  the  perturbation 
in  6H5  in  this  growing  mode  by  matching  at  (j)min  as  before.  One  finds  SH5  = 
12 xihf/c2  where  H§  is  the  background  value,  at  the  minimum.  Beyond  this 

point,  5H5  grows  as  f^/c  oc  eV^,  for  large  c,  whereas  in  the  background 
scaling  solution  decays  with  0  as  When  the  perturbation  is  of 

order  the  background  value,  the  trajectory  departs  from  the  scaling  solution 
and  the  potential  becomes  irrelevant.  The  departure  occurs  when  the  scalar 
held  has  attained  the  value 


0Dep 


2,  c2 
+  -In - , 

c  12*  ’ 


|V| 


< 


m 


V m 


(26) 


As  0  passes  beyond  0£>ep  the  kinetic  energy  overwhelms  the  negative  poten¬ 
tial  and  the  held  passes  onto  the  plateau  Vo  with  H5  nearly  constant  (see 
Figure  |2~2p,  and  equal  to 


s/6 


Ho(out)  ps  y 


12%, 


H5(in), 


(27) 


until  the  radiation,  matter  and  vacuum  energy  become  significant  and  i/5  is 
then  damped  away  to  zero. 

Before  moving  on  to  discuss  these  late  stages,  it  is  instructive  to  compare 
how  rapidly  0  travels  over  its  range  before  and  after  the  bounce.  The  time 
spent  to  the  left  of  the  potential  well  (0  <  0m,n)  is  essentially  identical  in 
the  incoming  and  outgoing  stages  for  x  <<  namely 
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For  the  outgoing  solution,  when  0  has  left  the  scaling  solution  but  before 
radiation  domination,  the  definition  Eq.  (16)  may  be  integrated  to  give  the 
time  since  the  big  bang  at  each  value  of  0, 


2  eVT^ 

3  H5(out ) 


(29) 


The  time  in  Eq.  (|29|)  is  a  microphysical  scale.  The  corresponding  formula  for 
the  time  before  the  big  crunch  is  very  different.  In  the  scaling  solution  (p^) 
one  has  for  large  c 


t(0) 


^min)/^ 
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min  )/2 
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(30) 


The  large  exponential  factor  makes  the  time  to  the  big  crunch  far  longer  than 
the  time  from  the  big  bang,  for  each  value  of  0.  This  effect  is  due  to  the 
increase  in  H§  after  the  bounce,  which,  in  turn,  is  due  to  the  positive  value 
of  X- 


5  The  Radiation,  Matter  and  Quintessence 
Epochs 

As  the  scalar  held  passes  beyond  the  potential  well,  it  runs  onto  the  positive 
plateau  V0.  As  mentioned  in  the  last  section,  the  value  of  H5(out )  is  nearly 
canceled  in  the  passage  across  the  potential  well,  and  is  reduced  to  H5  given 
in  Eq.  (|27|).  Once  radiation  domination  begins,  the  held  quickly  converges 
to  the  large  t  (Hubble-damped)  limit  of  Eq.  (pT|),  namely 

0c  =  ^ In  (2^H5(out)/Hr)^  ,  (31) 

where  Hr  is  the  Hubble  radius  at  kinetic-radiation  equality.  The  dependence 
is  obvious:  the  asymptotic  value  of  0  depends  on  the  ratio  of  H^out)  to 
Hr.  Increasing  H5(out)  pushes  0  further,  likewise  lowering  Hr  delays  radia¬ 
tion  domination  allowing  the  logarithmic  growth  of  0  in  the  kinetic  energy 
dominated  phase  to  continue  for  longer. 
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As  the  kinetic  energy  red  shifts  away,  the  gently  sloping  potential  gradu¬ 
ally  becomes  important,  in  acting  to  slow  and  ultimately  reverse  0’s  motion. 
The  solution  of  the  scalar  field  equation  is,  after  expanding  Eq.  ©  for  large 
r,  converting  to  proper  time  t  —  f  a (r) dr  and  matching, 


0 


\[?iHr 

a3(t) 


(32) 


where  as  above  we  define  a(t)  to  be  unity  at  kinetic-radiation  equal  density. 
During  the  radiation  and  matter  eras,  the  first  term  scales  as  t~  2  and  t~2  re¬ 
spectively.  For  a  slowly  varying  field,  is  nearly  constant,  and  the  potential 
gradient  term  in  Eq.  (^)  scales  linearly  with  t,  so  it  eventually  dominates. 

When  does  (j)  turn  around?  We  give  a  rough  discussion  here,  ignoring 
factors  of  order  unity.  First,  we  use  the  fact  that  V0  ~  t$  2  where  t0  is  the 
present  age  of  the  Universe,  and  roughly  we  have  to  ~  105tm,  where  tm  is 
the  time  of  matter  domination.  As  we  shall  see,  </>  may  reach  its  maximal 
value  (j)max  and  turn  around  during  the  radiation,  matter  or  quintessence 
dominated  epoch.  All  three  possibilities  are  acceptable  phenomenologically, 
although  the  case  where  turnaround  occurs  in  the  radiation  epoch  appears 
more  likely.  For  example,  (j)max  is  reached  in  the  radiation  era,  if,  from 
Eq.  fD, 


<  1. 


(33) 


If  the  Universe  becomes  radiation  dominated  at  the  GUT  scale,  tr  ~  1CU25 
seconds.  Then,  only  if  we  fine-tune  such  that  V/V'  >  1010  does  tmax  exceed 
tm.  This  corresponds  to  the  case  where  we  have  1010  e-foldings  or  more  of 
cosmic  acceleration  at  late  times  as  (f)  rolls  back,  far  more  than  required  for 
the  cyclic  solution.  The  bound  changes  somewhat  if  the  Universe  becomes 
radiation  dominated  as  late  as  nucleosynthesis  (tr  ~  1  second).  In  that 
case,  even  if  V /V^{{(f>c)  is  not  much  greater  than  unity,  the  scalar  field  turns 
around  in  the  matter  era  or  later.  For  turnaround  in  the  matter  era,  we 
require 

3  X  10-4  £  6  3  £  30.  (34) 


Finally,  if  the  field  runs  to  very  large  4>c,  so  that  V^/V(4>c)  ~  ce  C^c  is  expo¬ 
nentially  small,  then  (j)  only  turns  around  in  the  quintessence-dominated  era. 


For  the  example  considered  here,  the  natural  range  of  parameters  corresponds 
to  turnaround  occurring  during  the  radiation-dominated  epoch.  Hence,  by 
the  present  epoch,  the  field  is  rolling  monotonically  in  the  negative  direction 
and  slowly  gaining  in  speed.  Consequently,  the  ratio  of  the  pressure  to  the 
energy  density  is  increasing  from  its  value  at  turnaround,  w  —  —  1,  towards 
zero.  Depending  on  the  details  of  the  scalar  potential  V (0),  it  is  conceivable 
that  the  increasing  value  of  w  could  ultimately  be  observationally  detectable. 

Once  the  field  has  turned  around  and  started  to  roll  back  towards  the 
potential  well,  the  second  term  in  Eq.  fl32[)  dominates.  For  our  scenario 


to  be  viable,  we  require  there  to  be  a  substantial  epoch  of  vacuum  energy 
domination  (inflation)  before  the  next  big  crunch.  The  number  of  e-foldings 
Ne  of  inflation  is  given  by  the  usual  slow-roll  formula, 


Ne  = 


d(f> 


V_ 

VO 


?c(t>C 


(35) 


for  our  model  potential.  For  example,  if  we  demand  that  the  number  of 
baryons  per  Hubble  radius  be  diluted  to  below  unity  before  the  next  con¬ 
traction,  which  is  certainly  over-kill  in  guaranteeing  that  the  cyclic  solution 
is  an  attractor,  we  set  e3Ne  ^  1080,  or  Ne  ^  60.  This  is  easily  fulfilled  if  0c 
is  of  order  unity  in  Planck  units. 

From  the  formulae  given  above  we  can  also  calculate  the  maximal  value 
0c  dr  the  cyclic  solution:  for  large  c  and  for  tr  »  it  is 


0C  0mm  ~ 


(36) 


where  we  used  Hr  1  ~  tr,  the  beginning  of  the  radiation-dominated  epoch. 
From  Eq.  (^6|)  we  obtain 
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(37) 


This  equation  provides  a  lower  bound  on  tr.  The  extreme  case  is  to  take 
\Vmin\  ~  1.  Then  using  Vo  ~  10  120,  c  ~  10,  Ne  ~  60,  we  find  tr  ~ 
10~25  seconds.  In  this  case  the  maximum  temperature  of  the  Universe  is 
~  1010  GeV.  This  is  not  very  different  to  what  one  finds  in  simple  inflationary 
models. 
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As  (ft  rolls  down  the  hill,  one  can  check  that  (ft  leaves  the  slow-roll  regime 
when  e~c<^  exceeds  3/c2.  At  this  point  the  constant  term  Vo  in  the  potential 
becomes  irrelevant  and  one  can  use  the  scaling  solution  for  0,  all  the  way  to 
the  potential  minimum.  This  is  also  the  point  at  which  density  perturbations 
start  to  be  generated  via  the  ekpyrotic  mechanism,  while  the  Einstein  frame 
scale  factor  a  is  still  expanding.  The  Universe  continues  to  expand  slowly, 
but  with  a  slowly  decreasing  Hubble  constant,  and  finally  enters  contraction 
when  when  the  density  in  the  scalar  held  reaches  zero,  at  a  negative  value  of 
the  potential  energy.  The  ensuing  contracting  phase  is  accurately  described 
by  the  scaling  solution  (|22|),  in  which  a  ~  (—  t)p  and  <ft  =  y/2p/t,  with  t  <  0, 
and  t  =  0  being  the  time  of  the  next  big  crunch.  From  the  formulae  (|T^) 
one  finds  d\  =  ( pax  —  A Jp/3a0)/t ,  which  is  greater  than  zero  for  p  <  |,  since 
do  is  greater  than  a\.  Thus  even  when  a  is  undergoing  slow  contraction,  in 
the  scaling  era,  the  effect  of  the  motion  of  (ft  is  enough  to  make  cq  expand 
throughout  this  phase.  Matter  residing  on  this  brane  would  see  continuous 
expansion  all  the  way  to  the  big  crunch.  The  same  argument  shows  that  cq 
actually  undergoes  a  small  amount  of  contraction  in  the  very  much  shorter 
scaling  epoch  of  the  expanding  phase. 


6  Generation  of  Density  Perturbations 

In  the  cyclic  scenario,  the  period  of  exponential  expansion  occurring  late  in 

each  cycle  plays  a  key  role  in  diluting  the  densities  of  matter,  radiation  and 

black  holes  to  negligible  levels,  suppressing  long  wavelength  perturbations 

and  establishing  a  ‘clean  slate’,  namely  a  flat  vacuous  Universe  in  which  all 

fields  are  in  their  quantum  mechanical  ground  state.  As  the  scalar  held  rolls 

down  the  potential  in  Eq.  entering  the  scaling  solution  in  Eq.  (^2[),  the 

ekpyrotic  mechanism  for  the  generation  of  fluctuations  derived  in  Refs.  1 

and  [12]  sets  in  and  a  scale  invariant  spectrum  of  adiabatic  perturbations 

is  thereby  developed.  Quantum  fluctuations  of  the  usual  inhationary  sort 

are  also  developed  in  the  slow-roll  quintessence  epoch,  but  these  are:  (a) 

negligible  in  amplitude  because  Vo  is  tiny;  and,  (b)  only  excited  on  scales  of 

order  to  and  above  in  the  contracting  phase.  These  scales  are  shrunk  only  as 

(— t)p  in  the  contracting,  scaling  solution,  but  then  expanded  as  t a,  tz  and 
2 

1 3  in  the  kinetic  dominated,  radiation  and  matter  eras  in  the  big  bang  phase, 
which  also  lasts  for  a  time  of  order  to-  Therefore,  the  modes  amplified  during 
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inflation  are  exponentially  larger  in  wavelength  than  the  Hubble  radius  scale 
in  the  next  cycle  by  time  of  quintessence  domination,  which  is  the  present 
epoch. 

Let  us  concentrate  on  the  fluctuations  produced  via  the  ekpyrotic  mecha- 


the  inhomogeneous  fluctuations  in  the  scalar  field  S(j)(t,  x) 


msm.i 


Sg<50^(r)e*fc'‘T,  we  remove  the  damping  term  by  setting  50^  =  a  1x%,  to  obtain 


(38) 


where  primes  denote  conformal  time  derivatives  and  we  have  defined  kp, 
the  comoving  ‘freeze-out’  wavenumber.  Modes  with  k  >  kp  oscillate  with 
fixed  amplitude,  whereas  those  with  k  <  kp  are  amplified.  In  the  regime 
of  interest  kp  grows  monotonically  so  that  shorter  and  shorter  wavelengths 


progressively  freeze  out  as  the  big  crunch  is  approached.  The  physical  scale 
at  which  modes  freeze  out  is  given  by 


2 


(39) 


As  usual  we  adopt  units  where  8nG  =  1,  and  denote  proper  time  derivative 
with  dot.  In  the  era  of  quintessence  domination  when  V  dominates  over 
V^,  the  freeze-out  scale  Xp  is  nearly  constant,  and  comoving  wavelengths 
are  exponentially  stretched  beyond  it.  As  V ^  begins  to  dominate  however, 
Hubble  damping  becomes  irrelevant,  and  the  system  approaches  the  scaling 
solution  given  in  Eq.  (^2j),  in  which  ~  —2/t2,  where  t  is  the  proper 
time  to  the  big  crunch.  The  freeze-out  scale  drops  linearly  with  time  to 
zero,  as  the  scale  factor  is  falling,  like  (—t)s.  Therefore  progressively  shorter 
and  shorter  wavelength  modes  are  frozen  out  and  amplified,  with  waves  of 
physical  wavelength  tp  being  frozen  out  at  a  time  tp. 

An  exponentially  large  band  of  comoving  wavelengths  is  amplified  and 
frozen  in  as  0  rolls  from  0  =  0  down  towards  0mjn.  Modes  with  all  physical 
wavelengths  from  the  microphysical  scale  tmin ,  which  could  be  not  much 
larger  than  the  Planck  length,  to  the  macroscopic  scale  f0/c  which  is  of  order 
a  tenth  the  present  Hubble  radius,  aquire  scale  invariant  perturbations.  Once 
the  perturbations  are  generated,  their  wavelength  scales  as  (tmin/tp)v  in  the 
collapsing  phase.  Then  as  0  crosses  the  potential  well  and  races  off  to  minus 
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infinity,  the  Einstein  frame  physical  wavelength  goes  to  zero.  But  this  is  not 
the  relevant  quantity  to  track,  since  we  match  the  variables  ao  and  a\  and 
therefore  should  match  the  physical  wavelengths  as  measured  by  these  scale 
factors.  In  the  kinetic  dominated  phase,  a 0  and  a\  are  nearly  constant,  so  in 
effect  the  physical  wavelength  of  the  modes  are  matched  when  0  crosses  (ftmin , 
in  the  contracting  and  expanding  phases.  Furthermore,  the  contracting  and 
expanding  solutions  are  nearly  time-reverses  of  one  another,  until  the  time 
tjjep  computed  above  when  the  expanding  solution  deviates  from  scaling. 
Therefore  one  is  effectively  matching  at  toep,  from  which  one  sees  that  the 
time  ip  at  which  perturbations  on  the  current  Hubble  radius  scale  to  were 
generated,  is  given  by 
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(40) 


where  the  bracketed  factors  are:  (a)  the  contraction  of  the  scale  factor  in  the 
scaling  solution,  between  the  time  tp  at  which  the  perturbations  were  gener¬ 
ated  and  the  time  tpep  at  which  the  expanding  solution  departs  from  scaling; 
(b)  the  scaling  back  of  the  present  comoving  Hubble  radius  scale  to  the 
time  of  matter- domination  tm;  (c)  the  scaling  back  to  the  time  of  radiation- 
domination  tr]  and  (d)  the  scaling  back  to  the  time  toep  using  H5  ~  constant, 
corresponding  to  kinetic  domination  in  the  expanding  solution. 

From  Eq.  ([3(]),  it  follows  that  perturbations  on  the  scale  of  the  present 
Hubble  radius  were  generated  at  a  held  value 
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Comparing  with  Eq.  (^)  for  the  resting  value  of  the  held  0c,  and  the  ex¬ 
pression 
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for  the  held  value  at  the  potential  minimum,  which  follows  from  Eq.  (|28|), 
one  finds  that 
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where  the  first  term  dominates.  In  other  words,  the  fluctuations  we  see  today 
were  generated  at  a  held  value  approximately  half  way  between  the  zero  and 
the  minimum  of  U(0). 

7  Cyclic  Solutions  and  Cyclic  Attractors 

We  have  shown  that  a  cyclic  Universe  solution  exists  provided  we  are  allowed 
to  pass  through  the  Einstein-frame  singularity  according  to  the  matching  con¬ 
ditions  elaborated  in  Section  [|.  Eqs.  (^4j)  and  (^5]) .  Specifically,  we  assumed 
that  H5(out )  =  —  (1  +  x)H${in)  where  x  is  a  non- negative  constant,  cor¬ 
responding  to  branes  whose  relative  speed  after  collision  is  greater  than  or 
equal  to  the  relative  speed  before  collision.  Our  argument  showed  that,  for 
each  x  >  0,  there  is  a  unique  value  of  H5(out )  that  is  perfectly  cyclic.  In  the 
Appendix,  we  show  that  an  increase  in  velocity  is  perfectly  compatible  with 
energy  and  momentum  conservation  in  a  collision  between  a  positive  and 
negative  tension  brane,  provided  a  greater  density  of  radiation  is  generated 
on  the  negative  tension  brane.  (A  similar  outcome  can  occur  through  the 
coupling  of  0  to  the  matter  density,  as  discussed  below  Eq.  (|25|),  but  we  will 
only  discuss  the  Erst  effect  for  the  purpose  of  simplicity.) 

In  this  Section,  we  wish  to  show  that,  under  reasonable  assumptions,  the 
cyclic  solution  is  a  stable  attractor,  typically  with  a  large  basin  of  attraction. 
Without  the  attractor  property,  the  cyclic  model  would  seem  hne-tmied  and 
unstable.  One  could  imagine  that  there  would  still  be  brane  collisions  and 
periods  of  contraction  and  expansion,  but  there  would  be  no  regularity  or 
long-term  predictability  to  the  trajectories.  If  this  were  the  case,  fundamental 
physics  would  lose  its  power  to  explain  the  masses  and  couplings  of  elemen¬ 
tary  particles.  The  masses  and  couplings  depend  on  0  and  other  moduli 
holds.  If  there  were  no  attractor  solution,  the  precise  trajectory  of  0  through 
cosmic  history  would  depend  on  initial  conditions  and  could  not  be  derived 
from  fundamental  physics  alone.  In  our  proposal,  the  nature  of  the  attrac¬ 
tor  solution  depends  on  microphysics  at  the  bounce  which  is  computable,  in 
principle,  from  fundamental  theory.  Hence,  masses  and  couplings  of  particles 
change  during  the  course  of  cycle,  but  fundamental  theory  retains  predictive 
power  in  determining  the  way  they  change  and,  specifically,  their  values  at 
the  current  epoch. 

The  essential  feature  for  attractor  behavior  is  the  extended  period  of 
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accelerated  expansion  that  damps  the  motion  of  (j).  Let  us  consider  how 
this  works.  Assuming  y  is  fixed  by  microphysics,  there  is  a  value  H$(out) 
which  corresponds  to  the  cyclic  solution.  Now,  suppose  the  value  of  H5(out ) 
exceeds  H5(out).  This  means  that  the  outgoing  velocity  exceeds  the  cyclic 
value  and  (j)  runs  out  farther  on  the  plateau  than  in  the  cyclic  case.  Once  the 
field  stops,  turns  around,  and  quintessence-domination  begins,  the  field  is 
critically  damped.  By  the  time  V (</>)  falls  to  zero,  the  transient  behavior  of  </> 
which  depends  on  the  initial  value  of  H 5  has  damped  away  exponentially  so 
that  the  field  accurately  tracks  the  slow-roll  solution.  Following  the  solution 
forwards,  H^{in)  at  the  next  bounce  is  then  exponentially  close  to  what  it 
would  have  been  for  the  cyclic  solution.  By  erasing  memory  of  the  initial 
conditions,  the  acceleration  insures  that  H^out)  after  the  next  bounce  is 
very  nearly  H5(out). 

How  many  e-foldings  of  accelerated  expansion  are  actually  required  to 
make  the  cyclic  solution  an  attractor?  If  there  is  no  epoch  of  accelerated 
expansion,  perturbations  will  grow  each  cycle,  becoming  self-gravitating  and 
non-linear  so  that  no  attractor  will  occur.  A  minimal  requirement  for  ob¬ 
taining  an  attractor  is  that  linear  density  perturbations  grown  during  the 
matter  era  should  be  damped  away  during  the  subsequent  exponential  ex¬ 
pansion.  This  requires  at  least  ln(105)~  10  e-foldings  of  exponential  expan¬ 
sion.  Equally,  diluting  the  number  density  of  baryons  below  one  per  Hubble 
volume  is  certainly  over-kill  in  terms  of  ensuring  an  attractor,  and  this  re¬ 
quires  of  order  60  e-foldings.  In  fact,  as  we  discussed  above,  obtaining  a  far 
larger  number  of  e-foldings  is  perfectly  possible. 

To  discuss  the  nature  of  the  attractor  solution,  it  is  helpful  to  plot  the 
trajectories  of  the  system  in  the  phase  space  given  by  the  (if5,  0)-plane, 
shown  in  Figures  [|  through  |6|.  Recall  that  H5  is  proportional  to  </>;  see 
Eq.  (|1^).  Figure  [4]  illustrates  the  cyclic  trajectory  for  the  case  where  no 
radiation  is  generated  at  the  bounce  (y  =  0)  and  the  cycle  is  exactly  time- 
symmetric. 

The  phase  space  plot  must  always  satisfy  three  properties.  First,  for  a 
flat  Universe,  the  Friedmann  constraint  equation  H 2  =  \p  implies  that  the 
energy  density  p  =  (|</>2  +  U(</>))  must  be  positive.  Without  negative  space 
curvature,  the  system  is  simply  not  allowed  to  explore  negative  energies.  We 
show  the  classically  excluded  region  as  the  shaded  area  in  the  Figures.  In 
Figure  [|  where  there  is  no  radiation  and  U(0)  — >  0  as  0  — >  0,  the  excluded 
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region  extends  along  H5  —  0  out  to  0  —>  —  oo.  In  Figures  |5|  and  []  the 
shape  of  the  excluded  region  is  modified  due  to  the  presence  of  radiation. 
For  example,  the  grey  region  pinches  off  on  the  left  hand  side  for  some  finite 
value  of  0.  However,  the  effect  is  negligible  for  the  trajectories  considered  in 
our  discussion  and  so  we  show  the  same  excluded  region  in  the  Figures  as  in 
the  case  of  no  radiation. 

The  second  property  is  that  phase  space  trajectories  are  double-valued 
on  the  (i/5,  0)-plane.  Given  the  scalar  held  and  H^,  one  may  have  either  a 
contracting  or  an  expanding  Universe.  We  represent  expanding  trajectories 
as  solid  lines  and  the  contracting  trajectories  as  dashed  lines.  Two  expanding 
trajectories  are  not  allowed  to  cross,  and  neither  are  two  contracting  trajec¬ 
tories  for  the  usual  reasons  that  hold  for  particle  trajectories  on  phase  space. 
However,  an  expanding  trajectory  may  certainly  intersect  a  contracting  tra¬ 
jectory. 

The  final  rule  is  that  there  are  only  two  ways  an  expanding  trajectory  can 
turn  into  a  contracting  trajectory.  If  reversal  occurs  at  finite  0  it  can  only 
happen  if  the  trajectory  hits  the  forbidden  zero  density  region  (shaded),  since 
p  has  to  vanish  if  ii  is  to  pass  smoothly  through  zero.  The  shaded  region  is 
analogous  to  the  “egg”  region  described  by  Brustein  and  Veneziano^  The 
second  way  in  which  contraction  can  turn  into  expansion  is  if  the  system  runs 
off  to  0  =  —00.  Then,  the  ‘bounce’  described  in  Section  and  Ref.  |]  occurs. 

The  trajectory  shown  in  Figure  (|  is  a  cyclic  solution  (albeit  not  a  very 
interesting  one)  in  which  no  matter-radiation  is  produced  at  the  bounce  and 
the  value  of  H5(out)  is  precisely  equal  to  if5(m).  The  held  rolls  out  in 
the  expanding  phase  (solid  line  emanating  from  the  upper  left  side  of  the 
Figure)  to  the  value  where  U(0)  =0  and  stops  (the  rightmost  tip  of  the 
grey  region).  The  total  energy  density  is  momentarily  zero  and  expansion 
reverses  to  contraction.  The  held  then  rolls  back  to  —00  (lower  left  side  of 
the  Figure).  The  expanding  and  contracting  phases  are  exactly  symmetrical. 

The  time-scale  for  one  cycle  of  this  empty-Universe  solution  is  easily  esti¬ 
mated  by  noting  that  most  of  the  time  is  spent  near  the  zero  of  the  potential, 
where  it  is  rather  shallow,  and  the  scale  factor  a  is  nearly  constant.  There¬ 
fore,  we  can  neglect  gravity  and  calculate  the  period  for  one  cycle  in  the  case 
of  the  empty  Universe: 


-f-  r^~' 

0  empty 


1  rc\<t>min\  cty  1 

C\J\ 0  Jo  \Jev  —  1  Cy/V 0 


(44) 
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Figure  4:  The  cyclic  trajectory  in  the  (H $,(/))- plane  for  the  case  where 
no  matter  and  radiation  are  produced  at  the  bounce  (y  =  0).  The  grey 
region  which  corresponds  to  negative  energy  density,  is  forbidden.  The  solid 
(dashed)  line  represents  the  trajectory  during  an  expanding  (contracting) 
phase.  Expansion  turns  to  contraction  and  vice  versa  when  the  trajectory 
hits  the  zero  energy  surface  (the  rightmost  tip  of  grey  region  in  this  case). 
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for  large  c\(ftmin\.  For  the  parameters  typical  in  the  our  examples,  this  corre¬ 
sponds  to  roughly  one  tenth  of  the  current  age  of  the  Universe,  or  a  billion 
years. 

In  Figure  ||,  we  consider  the  case  where  radiation  is  produced  at  the 
bounce,  x  >  0-  If  H5(out)  is  too  low  compared  to  the  cyclic  value,  the 
trajectory  encounters  the  zero  density  boundary  (grey  region)  and  reverses 
to  contraction.  Solid  curves  represent  the  expanding  phase  of  the  trajectory, 
and  dashed  lines  represent  the  contracting  phase  of  the  trajectory.  Let  us 
call  this  an  “undershoot”  solution.  H 5  is  only  constant  if  (ft  and  a  are  both 
increasing,  or  both  decreasing.  However,  if  the  Universe  reverses  when  H§  is 
still  positive,  then  the  scalar  field  kinetic  energy  is  blue  shifted  and  and  H5 
is  rapidly  driven  to  more  and  more  positive  values.  The  trajectory  flies  off 
to  large  positive  H5  and  (ft  (the  upper  boundary). 

As  one  increases  H5(out),  the  behavior  of  the  system  changes.  For  suf¬ 
ficiently  large  H5l  the  system  avoids  the  zero  energy  surface  entirely  during 
the  period  when  (ft  is  increasing  (the  bold  solid  and  dashed  trajectory  in  Fig¬ 
ure  |5j) .  The  field  “overshoots”  the  negative  region  of  the  potential  and  lands 
on  the  positive  plateau.  Exponential  expansion  begins,  followed  by  a  very 
slow  roll  of  (ft  back  towards  the  potential  zero,  This  period  appears  as  a  long, 
thin  excursion  on  the  right  hand  side  of  the  figure.  \H^\  is  small  because  the 
field  is  rolling  slowing  in  the  quintessence-dominated  phase). 

The  cyclic  attractor  solution  lies  between  these  undershoot  and  overshoot 
trajectories.  Figure  shows  trajectories  with  initial  values  of  H5(out )  both 
above  and  below  the  cyclic  value  (the  middle  curve).  Here  we  can  study  the 
stability  of  the  cyclic  solution.  Let’s  first  consider  a  trajectory  with  H5(out) 
larger  than  the  value  in  the  cyclic  solution.  This  trajectory  is  indicated  by 
(1)  in  the  Figure.  Clearly,  it  overshoots  the  cyclic  trajectory  and  undergoes 
a  longer  period  of  exponential  expansion  (the  long  excursion  to  the  right). 
During  the  slow-roll  epoch,  the  difference  between  this  trajectory  and  the 
cyclic  one  damps  away  until  it  is  exponentially  small.  The  trajectory  en¬ 
counters  the  zero  density  surface  very  slightly  later  than  the  cyclic  solution 
does,  and,  therefore,  reverses  and  ends  up  with  a  very  slightly  smaller  value 
of  H^out)  than  that  in  the  cyclic  trajectory.  Similarly,  one  can  see  that 
starting  the  system  in  state  (2)  with  a  smaller  value  of  H5(out)  than  that  of 
the  cyclic  trajectory,  the  system  will  inflate  less  and  reverse  earlier,  ending 
up  with  a  larger  value  of  H5(out)  than  that  of  the  cyclic  trajectory.  This  dis¬ 
cussion  shows  that  the  trajectory  is  stable  and  that  memory  of  initial  decays 
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Figure  5:  Trajectories  in  the  (H5,  0)-plane  for  the  case  where  there  is  radia¬ 
tion.  The  solid  (dashed)  curves  represent  the  trajectory  during  an  expanding 
(contracting)  phase.  The  thin  lines  illustrate  undershoot  solutions  and  the 
heavy  line  represents  an  overshoot  solution. 
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exponentially  after  just  one  cycle. 

Here  we  implicitly  assumed  that  x  is  a  constant  independent  of  H^(in), 
the  incoming  velocity.  In  the  Appendix,  we  obtain  an  expression  for  x  hi 
Eq.  (^)  in  terms  of  the  matter-radiation  energy  densities  created  on  the 
positive  and  negative  tension  branes.  Assuming  the  energy  density  on  the 
negative  tension  brane  p_  is  significantly  greater  than  the  energy  density 
created  on  the  positive  tension  brane,  we  have  x  ^  P--  Assuming  the  col¬ 
lision  between  branes  occurs  at  a  low  velocity  so  that  one  is  not  far  from 
the  adiabatic  limit,  p_  should  decrease  with  decreasing  H5 in.  (Note  that  the 
low-velocity  assumption  has  been  made  throughout  since  it  is  required  for  the 
moduli  approximation.)  Hence,  we  would  anticipate  that  x  rises  monotoni- 
cally  as  the  incoming  velocity  increases.  This  effect  can  alter  the  trajectories 
and  the  precise  basin  of  attraction  in  detail,  but  does  not  alter  the  conclusion 
that  a  large  basin  of  attraction  exists.  This  is  assured  by  having  a  potential 
plateau  or,  more  generally,  a  region  of  the  potential  in  which  0  slow-rolls 
with  total  energy  comparable  to  the  current  dark  energy  density. 

Quantum  effects  are  also  unlikely  to  affect  the  attractor  solution.  We 
have  shown  that  solution  is  stable  under  small  perturbations  and  here  the 
perturbations  remain  small  since  |Emm|  and  Vo  are  small  compared  to  the 
Planck  scale. 

What  of  the  trajectories  in  Figure  |5|.  for  example,  that  run  away  to  large 
0?  For  these,  it  is  important  to  understand  what  happens  as  0  grows  more 
positive.  One  possibility  is  that  the  potential  F(0)  diverges  as  0  — >  oo.  Our 
example  for  F(0)  has  an  infinite  plateau,  but,  as  discussed  in  Section  |2j.  this 
is  not  a  general  requirement.  If  F(0)  grows  sufficiently,  0  will  bounce  back 
towards  —  oo.  Alternatively,  the  same  effect  can  occur  if  the  theory  includes 
massless  fields  that  couple  to  the  scale  factor  on  the  negative  tension  brane, 
a\.  The  Lagrangian  density  then  includes  a  term  oc  l/a\.  Increasingly 
positive  0  corresponds  to  shrinking  cq .  Hence,  this  dynamical  term  can  also 
create  a  force  that  causes  0  to  bounce  back.  The  net  effect  is  that  0  rattles 
back  and  forth  along  the  potential,  possibly  following  a  chaotic  orbit.H’H 
These  effects  could  enhance  the  basin  of  attraction  for  the  cyclic  solution. 
That  is,  some  of  these  trajectories  which  we  ignored  in  our  undershoot  and 
overshoot  treatment  may  eventually  hit  the  plateau  with  low  velocity,  at 
which  point  they  would  become  drawn  to  the  attractor  solution. 

Finally,  let  us  emphasize  that  we  have  only  considered  the  issue  of  sta¬ 
bility  in  the  context  of  the  very  simplified  model  studied  here,  with  a  single 
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Figure  6:  Trajectories  in  the  ( H 5,  0)-plane  showing  the  attractor  nature  of 
the  cyclic  solution  (the  middle  trajectory).  Path  (1)  is  an  overshoot  solution 
that  begins  with  slightly  greater  velocity  (H5)  than  the  attractor,  bounces 
off  the  grey  zero-energy  surface,  and  then  has  a  contracting  trajectory  whose 
value  of  H5  is  smaller  in  magnitude.  Path  (2)  is  an  undershoot  solution  which 
begins  with  slightly  less  velocity  in  the  expanding  phase  than  the  cyclic  value 
and  ends  in  a  contracting  phase  with  H§  having  a  slightly  greater  magnitude. 
Following  the  next  bounce,  therefore,  overshoot  turns  into  undershoot  and 
vice-versa.  In  either  case,  the  deviation  from  the  attractor  value  shrinks. 


40 


scalar  field  0,  and  the  matching  conditions  discussed  in  Section  [|.  It  would  be 
very  interesting  to  generalize  this  discussion  to  include  other  moduli,  matter 
which  couples  in  a  nontrivial  way  to  0  as  discussed  in  Section  and  and  also 
discrete  degrees  of  freedom  such  as  a  quantized  four-form  held,  which  may 
change  from  cycle  to  cycle  so  that  the  system  really  explores  moduli  space. 
The  existence  or  otherwise  of  an  attractor  could  well  be  relevant  to  the  deter¬ 
mination  of  the  relative  abundances  of  dark  energy,  dark  matter,  baryons  and 
photons  in  the  Universe,  and  also  to  the  values  of  the  fundamental  constants 
of  nature. 


8  Implications 

The  strengths  of  the  cyclic  model  are  its  simplicity,  its  efficient  use  of  all 
of  the  dominant  elements  of  the  Universe  and  the  fact  that  it  is  a  com¬ 
plete  description  of  all  phases  of  cosmic  evolution.  This  can  be  contrasted 
with  inhationary  cosmology,  a  highly  appealing  theoretical  model  in  its  own 
right.  Inhationary  cosmology  focuses  on  a  brief  epoch  when  the  Universe  was 
1CT35  seconds  old.  The  model  relies  on  assumptions  about  how  the  Universe 
emerged  from  the  cosmic  singularity.  One  must  postulate  the  existence  of 
a  phase  of  rapid  cosmic  acceleration  at  very  high  energies,  for  which  their 
exists  no  direct  proof.  (In  contrast,  the  cyclic  model  relies  on  low  energy  cos¬ 
mic  acceleration  that  has  been  observed.)  Subsequent  cosmic  events,  such 
as  the  recent  transition  from  matter-domination  to  dark  energy  domination 
and  cosmic  acceleration,  appear  to  have  no  direct  connection  to  inhationary 
theory. 

Because  the  cyclic  model  ties  the  past,  present  and  future  evolution  of  the 
Universe  in  a  tight,  cross-correlated  way,  it  has  surprising  explanatory  and 
predictive  power.  In  the  introduction  to  this  paper,  we  noted  a  number  of  the 
most  challenging  questions  of  cosmology  and  fundamental  physics.  In  this 
section,  we  consider  each  of  these  questions  (and  more)  and  briehy  describe 
the  insights  the  cyclic  model  provides  concerning  their  answers. 
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8.1  Why  is  the  Universe  homogeneous,  isotropic  and 
flat? 

The  Universe  is  made  homogeneous  and  isotropic  during  the  period  of  the 
preceding  cycle  when  quintessence  dominates  and  the  Universe  is  undergoing 
slow  cosmic  acceleration.  This  ensures  that  the  branes  are  flat  and  parallel 
as  they  begin  to  approach,  collide,  and  emerge  in  a  big  bang.  Inflation  also 
relies  on  cosmic  acceleration,  but  driven  by  very  high  vacuum  energy  which 
produces  an  acceleration  that  is  nearly  lO100  times  faster. 

8.2  How  were  density  inhomogeneities  generated? 

In  the  cyclic  model,  the  observed  inhomogeneities  in  the  Universe  are  gener¬ 
ated  during  the  contracting  phase  when  the  scale  factor  is  nearly  static  and 
gravitational  effects  are  weak.  Consequently,  as  in  the  ekpyrotic  scenario,  a 
nearly  scale  invariant  spectrum  of  adiabatic,  gaussian  energy  density  fluctu¬ 
ations  is  generated.  However,  because  the  expansion  rate  is  negligible  and 
gravitational  effects  are  weak,  the  tensor  (metric  fluctuation)  spectrum  is 
blue  with  an  exponentially  tiny  amplitude  at  long  wavelengths. 

Fluctuations  are  also  created  during  the  quintessence  dominated  phase, 
just  as  they  are  during  inflation.  However,  because  the  energy  density  during 
the  accelerating  phase  is  100  orders  of  magnitude  smaller  than  in  inflation,  the 
resulting  fluctuation  amplitude  is  exponentially  smaller  in  the  cyclic  model. 
These  fluctuations  also  have  wavelengths  that  exceed  the  current  Hubble 
horizon.  Hence,  they  are  observationally  irrelevant. 

8.3  What  is  the  role  of  dark  energy  and  the  current 
cosmic  acceleration? 

Clearly,  dark  energy  and  the  current  cosmic  acceleration  play  an  essential 
role  in  the  cyclic  model  both  by  reducing  the  entropy  and  black  hole  density 
of  the  previous  cycle,  and  triggering  the  turnaround  from  an  expanding  to 
a  contracting  phase.  (In  all  other  cosmologies  to  date,  including  inflation, 
dark  energy  has  no  essential  role.) 
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8.4  How  old  is  the  Universe? 


A  truly  cyclic  Universe  is  clearly  infinitely  old  in  terms  of  cosmic  time.  As 
we  have  noted,  the  exact  cyclic  solution  can  also  be  an  attractor.  Hence, 
the  cycling  is  stable.  Consequently,  one  becomes  insensitive  to  the  initial 
conditions  for  the  Universe  as  long  as  they  were  within  the  basin  of  attraction 
of  the  cyclic  solution.  We  believe  that  within  this  framework,  the  problem  of 
the  initial  conditions  for  the  universe  is  significantly  altered:  as  long  as  the 
universe  has  some  nonzero  probability  for  entering  the  cyclic  solution,  large 
regions  of  the  universe  maintain  cyclic  evolution  for  arbitrarily  long  periods 
of  time. 

There  is  a  possible  objection  to  this  argument,  due  to  the  fact  that  the  four 
dimensional  nonsingular  brane  spacetimes  in  our  scenario  are  past  geodesi- 
cally  incomplete.  As  we  have  explained,  for  most  of  cosmic  time  they  are 
well  approximated  by  de  Sitter  spacetime,  with  a  cosmological  constant  (or 
vacuum  energy)  close  to  the  currently  observed  value.  This  nearly  de  Sitter 
spacetime  is  foliated  by  slices  of  constant  scalar  held  which  are  nearly  ge¬ 
ometrically  hat.  Matter  is  repeatedly  generated  on  the  slices  with  <f>  =  — oo, 
in  the  rest  frame  defined  by  those  slices. 

As  one  follows  cosmic  time  t  backwards,  one  must  pass  an  inhnite  num¬ 
ber  of  these  big  crunch-big  bang  surfaces.  However,  even  though  the  cosmic 
time  tends  to  —  oo,  the  proper  time  as  measured  along  timelike  geodesics 
running  into  the  past  generically  is  finite  even  as  t  tends  to  —  oo.  This  may 
be  seen  as  follows.  Consider  a  particle  with  momentum  P  in  the  hat  slicing. 
Its  momentum  blueshifts  as  a-1  as  you  follow  it  back  in  time.  The  geodesic 
becomes  nearly  null  and  the  proper  time  measured  along  the  geodesic  con¬ 
verges  even  though  t  tends  to  —  oo  (this  is  the  crux  of  the  recent  argument  of 
Borde,  Guth  and  Vilenkin  that  inflation  is  past  geodesically  incompleteEl) . 

In  our  scenario,  however,  all  physical  particles  are  created  with  finite  mo¬ 
mentum  in  the  flat  slicing  defined  by  fi.  If  we  follow  a  particle  present  in 
today’s  universe  back  in  time,  most  likely  it  was  created  on  the  last  (j)  =  —  oo 
surface.  With  an  exponentially  smaller  probability,  it  could  have  been  cre¬ 
ated  on  the  penultimate  =  —  oo  surface,  and  so  on  into  the  past.  The 
probability  that  any  observed  particle  originated  on  the  t  =  —  oo  hat  sur¬ 
face  which  is  the  boundary  of  the  hat  slicing  of  de  Sitter  spacetime,  is  zero. 
Therefore  we  do  not  attribute  any  physical  significance  to  the  past  geodesic 
incompleteness  of  the  spacetime  metric  in  our  scenario.  In  contrast,  particle 
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production  in  standard  inflationary  models  occurs  on  open  slices  a  finite  time 
ago. 

Even  if  there  are  no  particles  present  which  ‘saw’  the  past  boundary  of 
the  cyclic  universe,  one  might  object  that  the  scenario  implicitly  requires 
a  boundary  condition  in  the  infinite  past.  We  do  not  think  this  is  a  strong 
objection.  If  the  cyclic  solution  were  begun  within  a  finite  region  (for  example 
a  torus)  of  three  dimensional  space,  it  would  grow  exponentially  with  each 
cycle  to  an  arbitrarily  large  size.  After  an  arbitrarily  long  time,  to  any  real 
observer  the  universe  would  appear  to  be  infinite  both  in  spatial  extent,  and 
in  lifetime  to  the  past. 

So,  whilst  the  cyclic  model  still  requires  an  initial  condition,  provided  that 
state  is  within  the  basin  of  attraction  of  the  cyclic  solution,  we  are  completely 
insensitive  to  its  details.  Any  features  of  the  initial  state  (the  total  size  of 
the  Universe,  or  any  fluctuations  about  flatness  or  homogeneity),  become 
exponentially  diluted  in  each  cycle  and  since  the  cyclic  solution  can  repeat 
forever,  are  ultimately  completely  irrelevant  to  any  observation. 

8.5  What  is  the  ultimate  fate  of  the  Universe? 

The  cycles  can  be  continued  to  the  infinite  future,  as  well  as  the  infinite  past. 
Hence,  the  Universe  endures  forever. 

8.6  How  big  is  the  Universe? 

From  the  effective  4d  point-of-view,  the  Universe  oscillates  between  periods 
of  expansion  to  periods  of  contraction  down  to  a  big  crunch.  However,  from 
the  brane  world  point-of-view,  the  Universe  is  always  infinite  in  the  sense  that 
the  branes  always  have  infinite  extent.  The  fact  that  the  branes  are  spatially 
infinite  means  that  it  is  possibly  for  the  total  entropy  in  the  Universe  to 
increase  from  cycle  to  cycle,  and,  at  the  same  time,  have  the  entropy  density 
(in  particular,  the  total  entropy  per  Hubble  horizon)  become  nearly  zero  prior 
to  each  bounce. 

8.7  What  occurs  at  the  big  bang  singularity? 

The  cyclic  model  utilizes  the  ekpyrotic  notion  that  the  singularity  corre¬ 
sponds  to  the  collision  and  bounce  of  two  outer  orbifold  branes  in  a  manner 
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that  is  continuous  and  well-behaved.  The  singularity  is  not  a  place  where 
energy  and  curvature  diverge  and  time  begins.  Rather,  formulated  in  ap¬ 
propriate  fields  and  coordinates,  the  singularity  is  a  smooth,  finite  transition 
from  a  contracting  phase  heading  towards  a  big  crunch  and  a  big  bang  evolv¬ 
ing  into  an  expanding  Universe. 

8.8  What  determines  the  arrow  of  time? 

Since  the  Universe  is  cyclic,  it  may  appear  that  there  is  no  well-defined  means 
of  determining  the  arrow  of  time.  Indeed,  for  a  local  observer,  there  is  no 
clear  means  of  doing  so. 

From  the  global  perspective,  though,  there  is  a  clear  means  of  determining 
forward  from  backward  in  time.  First,  one  of  the  boundary  branes  is  forever 
expanding  in  the  “forward”  time  direction  in  the  cyclic  model.  The  other 
brane  is  expanding  except  for  brief  intervals  of  contraction,  but,  averaged 
over  a  cycle,  the  net  effect  is  expansion.  The  rate  changes  from  phase  to 
phase,  as  well  as  the  separation.  In  the  contraction  phase,  the  branes  them¬ 
selves  stretch  at  a  rate  that  is  slow  but  their  separation  rapidly  decreases.  In 
the  radiation,  matter,  and  quintessence  dominated  phases,  the  branes  stretch 
significantly,  but  their  separation  remains  fixed.  During  this  period,  the  en¬ 
tropy  created  during  the  previous  cycle  is  spread  out  exponentially,  reducing 
the  degrees  of  freedom  per  horizon  to  nearly  zero. 

8.9  Why  is  the  cosmological  constant  so  small? 

The  cyclic  model  provides  a  fascinating  new  outlook  on  this  vexing  problem. 
Historically,  the  problem  is  assumed  to  mean  that  one  must  explain  why  the 
vacuum  energy  of  the  ground  state  is  zero. 

In  the  cyclic  model,  the  vacuum  energy  of  the  ground  state  is  not  zero. 
It  is  negative  and  its  magnitude  is  large,  as  is  obvious  from  Figure  1.  If 
the  Universe  begins  in  the  ground  state,  the  negative  cosmological  constant 
will  cause  rapid  recollapse,  as  expected  for  an  anti-de  Sitter  phase.  In  the 
cyclic  scenario,  though,  we  have  shown  how  to  arrange  conditions  where  the 
Universe  avoids  the  ground  state.  Instead,  the  Universe  hovers  from  cycle  to 
cycle  above  the  ground  state  bouncing  from  one  side  of  the  potential  well  to 
the  other  but  spending  most  time  on  the  positive  energy  side.  The  branes 
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are  moving  too  rapidly  whenever  the  separation  corresponds  to  the  potential 
minimum. 

There  remains  the  important  challenge  of  explaining  why  the  the  cur¬ 
rent  potential  energy  is  so  small.  The  value  depends  on  both  the  shape  of 
the  potential  curve  and  the  precise  transfer  of  energy  and  momentum  at  the 
bounce.  Perhaps  explaining  the  value  will  be  an  issue  as  knotty  as  the  cos¬ 
mological  constant  problem,  or  perhaps  the  conditions  will  prove  easier  to 
satisfy.  What  is  certain,  though,  is  that  the  problem  is  shifted  from  tuning  a 
vacuum  energy,  and  this  provides  an  opportunity  for  new  kinds  of  solutions. 

8.10  Equation-of-state  of  dark  energy 

The  equation-of-state  of  the  dark  energy,  w,  is  the  ratio  of  the  pressure  to 
the  energy  density  of  0,  (|02  —  H)/(|02  +  V ).  In  Section  j|,  we  discussed  the 
evolution  of  0  in  the  radiation,  matter  and  quintessence  dominated  epochs. 
The  generic  result  is  that  evolution  of  0  in  the  positive  direction  halts  and 
the  field  begins  to  roll  back  towards  — oo  in  the  radiation-dominated  epoch. 
At  the  turn-around,  w  =  —  1  since  the  kinetic  energy  is  zero.  As  the  field  rolls 
back  and  its  kinetic  energy  increases,  w  increases.  Hence,  the  generic  result 
is  that  w  is  close  to  -1  today  and  increasing.  Conceivably,  cosmological 
observations  could  detect  this  prediction.  Tracker  models  of  quintessence, 
some  of  the  best-motivated  alternatives,  have  the  opposite  trend:  w  is  near 
—0.8  or  so  today  and  decreasing  towards  —1.0  Other  models,  such  as  k- 
essence,  have  the  same  trend  as  found  in  the  cyclic  mo  del. 0 


8.11  Implications  for  Supersymmetry  and  Superstrings 

The  cyclic  model  imposes  different  constraints  on  fundamental  physics  com¬ 
pared  to  previous  cosmological  models.  As  an  example,  consider  the  problem 
of  designing  supergravity  potentials.  The  potentials  are  constructed  from  a 
superpotential  W  according  to  the  prescription: 


-  „K/Ml 


V  =  e 


KijDiW  DAY 


M2pl 


WW 


(45) 


where  D,  =  d/dcjA  +  Ki/M^  is  the  Kahler  covariant  derivative,  Kt  =  dK/dcj)1, 
K.ij  =  d2K/ and  a  sum  over  each  superfield  0*  is  implicit.  If  the  ground 
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state  is  supersymmetric,  DiW  =  0,  the  first  term  is  zero.  In  general,  unless 
W  is  zero  for  precisely  the  same  values  for  which  DiW  =  0,  the  minimum 
has  a  negative  cosmological  constant.  In  the  past,  this  type  of  model  would 
have  been  ruled  unacceptable.  The  possibility  of  a  cosmology  in  which  the 
Universe  hovers  over  the  ground  state  in  a  state  of  zero  or  positive  energy 
revives  these  models  and  alters  constraints  on  model  building. 

An  obvious  but  important  implication  is  that  supersymmetry  breaking 
can  be  achieved  without  having  spontaneous  symmetry  breaking  in  the  ground 
state.  In  this  scenario,  it  suffices  if  the  Universe  hovers  in  the  radiation,  mat¬ 
ter  and  quintessence  dominated  epochs  at  some  state  far  above  the  ground 
state  in  energy  and  that  the  supersymmetry  is  broken  by  the  appropriate 
amount  in  the  hovering  state,  where  the  radiation,  matter  and  quintessence 
dominated  phases  occur.  These  considerations  have  a  significant  impact  on 
the  design  of  phenomenological  supersymmetric  models. 

One  other  requirement/prediction  of  the  cyclic  scenario  (and  the  ekpy- 
rotic  models  in  general)  is  that  the  branes  move  in  a  space-time  with  co- 
dimension  one.  The  constraint  derives  from  having  a  bounce  that  produces 
a  smooth  transition  from  contraction  to  expansion.  As  argued  by  Khoury 
et.  al,  the  geometry  is  flat  arbitrarily  close  to  the  bounce  provided  there  is 
one  extra  dimension  only.  Hence,  brane  world  scenarios  based  on  theories 
like  that  of  Horava  and  Witten  are  acceptable,  but  large  extra-dimensional 
models  relying  on  having  co-dimension  two  or  greater  are  problematic. 

8.12  Hoyle’s  Revenge? 

Within  each  cycle,  there  is  a  sequence  of  kinetic  energy,  radiation,  matter  and 
quintessence  dominated  phases  of  evolution  that  are  in  accord  with  the  stan¬ 
dard  big  bang  cosmology.  However,  averaged  over  many  cycles,  the  model  can 
be  viewed  as  a  remarkable  re-incarnation  of  Fred  Hoyle’s  steady  state  model 
of  the  Universe.  Most  of  the  cycle  is  spent  in  a  phase  with  nearly  constant 
energy  density,  as  in  the  steady  state  picture.  Indeed  Hoyle’s  C-held  that 
was  introduced  to  provide  a  constant  supply  of  matter  (and  a  preferred  rest 
frame)  is  replaced  by  our  scalar  held  0,  which  defines  a  preferred  time  slicing 
and  generates  matter  repeatedly  at  each  bounce,  restoring  the  Universe  to  a 
state  of  high  temperature  and  matter  density.  In  Hoyle’s  steady  state  model, 
every  hat  spatial  slice  was  statistically  identical.  Here  the  slices  are  iden¬ 
tical  only  when  separated  by  one  period,  so  we  have  a  discrete  rather  than 
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continuous  time  translation  symmetry.  Nevertheless  when  coarse  grained 
over  large  time  spans,  the  structure  is  similar  to  that  proposed  in  the  steady 
state  Universe.  Global  properties  of  the  cyclic  cosmology  will  be  discussed 
elsewhere  .0 
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9  Appendix:  Matching  H-}  across  the  Bounce 

In  this  appendix  we  discuss  the  matching  condition  needed  to  determine 
Hrj(out)  in  terms  of  H^{in).  We  shall  assume  that  all  other  extra  dimensions 
and  moduli  are  fixed,  and  the  bulk  space-time  between  the  branes  settles 
down  to  a  static  state  after  the  collision.  (In  the  simplest  brane  world  models, 
there  is  a  Birkhoff  theorem  which  ensures  that  there  is  a  coordinate  system 
in  which  the  bulk  metric  is  static  in  between  the  branes).  We  shall  take  the 
densities  of  radiation  on  the  branes  after  collision  as  being  given.  By  imposing 
Israel  matching  in  both  initial  and  final  states,  as  well  as  conservation  of  total 
energy  and  momentum,  we  shall  be  able  to  completely  fix  the  state  of  the 
outgoing  branes  and  in  particular  the  expansion  rate  of  the  extra  dimension 
Hs(out),  in  terms  of  H§{in).  A  more  complete  discussion  of  this  method  will 
be  presented  in  Ref. 

The  idea  is  to  treat  the  brane  collision  as  a  short-distance  phenomenon. 
The  warp  factor  may  be  treated  as  linear  between  the  branes  as  they  approach 
or  recede.  Linearity  plus  Zj  symmetry  ensures  that  the  kinks  in  the  warp 
factors  are  equal  in  magnitude  and  opposite  in  sign.  Israel  matching  relates 
the  kink  magnitudes  to  the  densities  and  speeds  of  the  branes,  yielding  the 
relations  we  use  below. 

The  initial  state  of  empty  branes  with  tensions  T  and  — T,  and  with 
corresponding  velocities  v+  <  0  and  V-  >  0  (measured  in  the  frame  in  which 
the  bulk  is  static)  obeys 

T\/l  -  v%  =  T/l  - 
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T 


Etot 

Ptot 


Tv+  Tv _ 

\A  !,+  V1  - 


(46) 


The  first  equation  follows  from  Israel  matching  on  the  two  branes  as  the 
approach,  and  equating  the  kinks  in  the  brane  scale  factors.  The  second  and 
third  equations  are  the  definitions  of  the  total  energy  and  momentum.  The 
three  equations  (f|6|)  imply  that  the  incoming,  empty  state  has  v+  =  — u_, 
Etot  =  0  and  that  the  total  momentum  is 


Ptot  — 


T  LHP  in) 


i/iqww 


<0, 


(47) 


where  we  identify  v+  —  V-  with  the  contraction  speed  of  the  fifth  dimension, 
\LHs(in)\. 

The  corresponding  equations  for  the  outgoing  state  are  easily  obtained, 
by  replacing  T  with  T  +  p+  =  T+  for  the  positive  tension  brane,  and  —  T 
with  —  T  +  p_  =  — T_  for  the  negative  tension  brane,  assuming  the  densities 
of  radiation  produced  at  the  collision  on  each  brane,  p+  and  p_  respectively, 
are  given  from  a  microphysical  calculation,  and  are  both  positive. 

We  now  wish  to  apply  energy  and  momentum  conservation,  and  Israel 
matching  to  the  final  state.  The  only  subtlety  is  that  the  (t,  y)  frame  in 
which  the  bulk  is  static  is  not  necessarily  the  same  frame  in  the  final  state 
as  it  was  in  the  initial  state,  so  one  should  boost  the  initial  two-momentum 
(Etot,  Ptot )  with  a  velocity  V  and  then  apply  the  Israel  constraints  and  energy- 
momentum  conservation  equations  in  the  new  boosted  frame.  The  latter 
provide  three  equations  for  the  three  unknowns  in  the  final  state,  namely 
v+(out),  V-(out)  and  V.  Writing  v± -(out)  =  tanh(0±),  where  Q±  are  the 
associated  rapidities,  one  obtains  two  solutions 
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where  T+  =  T+p+,  TL  =  T  —  p_  with  p+  and  p_  the  densities  of  radiation  on 
the  positive  and  negative  tension  branes  respectively,  after  collision.  Both  p+ 
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and  p-  are  assumed  to  be  positive.  In  the  first  solution,  with  signs  (— +),  the 
velocities  of  the  positive  and  negative  tension  branes  are  the  same  after  the 

collision  as  they  were  before  it.  In  the  second,  with  signs  ( - ),  the  positive 

tension  brane  continues  in  the  negative  y  direction  but  the  negative  tension 
brane  is  also  moving  in  the  negative  y  direction. 

The  corresponding  values  for  v±(out )  and  V  are 


v+(out ) 


V-  (out) 


V 


or 


\Ptot\  +  \Ptot\ 
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where  the  first  solution  for  V  holds  for  the  ( — h)  case,  and  the  second  for  the 
( - )  case. 

We  are  interested  in  the  relative  speed  of  the  branes  in  the  outgoing 
state,  since  that  gives  the  expansion  rate  of  the  extra  dimension,  —  v+(out)  + 
V-(out)  =  LH5(out),  compared  to  their  relative  speed  —  2v+  =  —LH5(in)  in 
the  incoming  state.  We  find  in  the  ( — h)  solution, 


H5(out) 

v+(out)  —  V-(out) 

Pit  +  4  T2 

H5(in) 

2v+  \ 

Pit  +  2 (T2  +  Tl)  +  Pt~J(Tl  -  Tl)2 

(50) 


and  in  the  ( - )  solution 


H5(out ) 
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Pit  +  4T2 

H5(in) 
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Pit  +  2 (T2  +  Tl)  +  Pt~2(T2  -  Tl)2 

with  Ptot  given  by  ([47])  in  both  cases. 

At  this  point  we  need  to  consider  how  the  densities  of  radiation  p+  and  p_ 
depend  on  the  relative  speed  of  approach  of  the  branes.  At  very  low  speeds, 
\LH5(in)\  «  1,  one  expects  the  outer  brane  collision  to  be  nearly  adiabatic 
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and  an  exponentially  small  amount  of  radiation  to  be  produced.  The  (— +) 
solution  has  the  speeds  of  both  branes  nearly  equal  before  and  after  collision: 

we  assume  that  it  is  this  solution,  rather  than  the  ( - )  solution  which  is 

realised  in  this  low  velocity  limit. 

As  \LH5(in)\  is  increased,  we  expect  p+  and  p_  to  grow.  Now,  if  we 
consider  p+  and  p_  to  be  both  <<  Ptot  «  T,  then  the  second  term  in 
the  denominator  dominates.  If  more  radiation  is  produced  on  the  negative 
tension  brane,  p_  >  p+,  then  \H5(out) / H5(in)\  =  (1  +  y)  ~  (1  +  (p~  — 
p_i_)/2T)  and  so  x  is  small  and  positive.  This  is  the  condition  noted  in  the 
text,  necessary  to  obtain  cyclic  behavior.  Conceivably,  the  brane  tension 
can  change  from  T  to  T'  =  T  —  t  at  collision.  Then,  we  obtain  (1  +  x)  ~ 
(1  +  (p_  —  p+  +  2t)/2T) 

For  the  ( — b)  solution,  we  can  straightforwardly  determine  an  upper  limit 
for  \H5(out)  /  H5(in)\  =  (1  +  y).  Consider,  for  example,  the  case  there  the 
brane  tension  is  unchanged  at  collision,  t  =  0.  The  expression  in  O  gives 
\H5(out)/H5(in)\  as  a  function  of  T+,  T!_  and  Ptot ■  It  is  greatest,  at  fixed  T_ 
and  Ptot. ,  when  T+  =  T,  its  smallest  value.  For  P2ot  <  T2,  it  is  maximized 
for  T2  =  T2  —  P2ot,  and  equal  to  \J  1  +  P2ot/ (4 T2)  when  equality  holds.  For 
Ptot  >  T2,  it  is  maximized  when  TL  =  0,  its  smallest  value,  and  P2ot  =  2T 2, 
when  it  is  equal  to  This  is  more  than  enough  for  us  to  obtain  the  small 
values  of  y  needed  to  make  the  cyclic  scenario  work.  A  reduction  in  brane 
tension  at  collions  t  >  0  further  increases  the  maximal  value  of  the  ratio. 
To  obtain  cyclic  behavior,  we  need  y  to  be  constant  from  bounce  to  bounce. 
That  is,  compared  to  the  tension  before  collision,  the  fractional  change  in 
tension  and  the  fractional  production  of  radiation  must  be  constant. 

We  shall  not  consider  the  ( - )  solution  in  detail,  except  to  note  that  in 

the  small  Ptot  limit  it  allows  an  arbitrarily  large  value  for  \H^{out) / H^{in)\, 
which  seems  unphysical. 

Let  us  re-iterate  that,  there  are  many  caveats  attached  to  this  calculation. 
We  have  not  calculated  p+  and  p_  and  have  left  these  as  parameters.  We 
have  set  the  brane  tensions  T  to  be  equal  before  and  after  the  collision, 
and  have  neglected  possible  bulk  excitations,  treating  the  bulk  spacetime  as 
static  both  before  and  after  collision.  We  have  ignored  matter  couplings  to 
bulk  scalars,  and  ignored  the  possible  dynamical  evolution  of  additional  extra 
dimensions.  Nevertheless  we  think  it  encouraging  that  the  unusual  behavior 
of  matter  bound  to  a  negative  tension  brane  allows  H5(out) / H5(in)  to  be 
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slightly  greater  than  unity,  which  is  what  we  need  for  cyclic  behavior. 
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